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ABSTRACT 


A  design  concept  for  a  non-rigid  energy -storage  flywheel  suitable  for 
home  or  farm  use  has  been  investigated  at  Texas  A§M  University  under  a  Sandia 
Laboratories'  contract.  Tlie  distinguishing  featiare  of  this  "flexible 
fl>nArheel"  is  its  construction  from  high  strength  fibers  (such  as  synthetic 
rope)  with  no  bonding  agent.  The  work  done  to  date  indicates  that  the 
flexible  flyivheel  is  self-balancing,  safe,  has  a  high  energy -density  capabil¬ 
ity  (60  Wh/lb) ,  and  should  be  simple  and  economical  to  manufacture.  The 
major  thrust  of  the  contract  work  has  been  to  find  a  solution  to  the 
subsynchronous  whirling  instability,  a  problem  common  to  all  high  speed 
flywheels,  which  is  especially  severe  in  the  flexible  flywheel  due  to  the 
high  supercritical  speed  ratio  and  the  large  amount  of  internal  friction. 

A  unique  gimbal  support  system  has  been  designed  and  analyzed  and  tested  on 
this  contract  which  stabilizes  the  flywheel  without  the  need  for  a  squeeze 
film  damper.  The  conceptual  design  was  developed  for  a  flexible  flywheel 
energy  storage  system  suitable  for  interfacing  with  a  small-scale  solar 
energy  source.  Cost  estimates  were  prepared  for  the  system  in  the  10  KWh 
and  50  KWh  sizes. 
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INTRODUCTION 


In  1975,  Dr.  R.  T,  Schneider  at  the  University  of  Florida  conceived  the 

idea  of  a  flexible  flyviieel  made  of  rope  for  energy  storage.  The  idea  was 
to  develop  an  economical,  safe,  self-balancing  energy  storage  device  to 
make  solar  or  wind -generated  electricity  practical  for  home  or  farm  use . 

Beginning  in  1976,  a  flywheel  test  facility  was  constructed  and  the 
first  rope  flyw'heels  were  spun  up.  It  soon  became  evident  that  one  of  the 
advantages  of  the  flexible  flywheel,  its  self-balancing  feature,  had  been 
bought  at  the  price  of  a  subs\Tichronous  whirling  instability  caused  by 
internal  friction.  This  poses  a  problem  because  the  whirl  critical  speed 
is  well  below  the  operating  speed  range  for  the  fl)nvheel.  Tlie  early  test 
work  on  the  flexible  fly\\iieel  was  carried  out  in  the  Nuclear  Engineering 
Department  at  the  University  of  Florida  by  Dr.  Schneider  and  several  other 
nuclear  engineering  students  who  had  little  or  no  experience  or  training  in 
the  field  of  rotor  dynamics.  In  fact,  it  was  not  realized  at  first  that 
the  whirling  was  not  synchronous  and  therefore  not  affected  by  rotor 
unbalance.  One  of  the  early  unsuccessful  efforts  to  control  the  whirling 
was  by  increasing  the  support  rope  constraints,  an  approach  which  rotor 
dynamics  theory  shows  will  actually  aggravate  the  problem. 

Since  the  major  tecimical  problem  of  the  flexible  flywheel  was  found 
to  be  one  of  rotor  dynamics,  the  author  of  this  report  (who  is  active  in 
this  field)  joined  work  on  the  project  while  in  the  Department  of  Mechanical 
Engineering  at  the  University  of  Florida.  Tlie  project  was  subsequently 
moved  to  Texas  A§M  University  in  1978  v/hen  the  author  moved  there,  and  all 
exjDeriments  and  analyses  described  in  this  report  were  performed  there . 
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'file  flexible  flywheel  concept  originated  by  Dr.  R.  T.  Schneider  at 
the  University  of  Florida  has  been  investigated  experimentally  and  ana- 
lyti.cally  at  Texas  Am  University.  Tlie  flexible  flywheel  is  constructed 
From  swthetic  rope  with  no  bonding  agent,  and  is  flexibly  supported  (no 
hub)  to  induce  self-balancing  characteristics. 

Tliree  experimental  flywheels  were  built  and  tested;  a  10-1/2"  dia¬ 
meter  N'ylon  flwheel  weighing  2.7  lb.,  a  22-1/2"  diameter  Dacron  wheel 
weighing  19  lb.,  and  a  30"  diameter  Nylon  wheel  weighing  88  lb.  All 
three  wheels  were  suspended  by  various  support  rope  (and/or  wire  cable) 
configurations  from  a  gimballed  motor/ generator  support  system  designed 
to  suppress  subsynchronous  whirl.  In  all  cases  the  motor/generator  was 
operated  in  the  motor  mode  only. 

The  experiments  demonstrated  that  the  gimballed  motor/generator 
support  system  has  a  strong  stabilizing  influence  on  flexible  flywheel  dy¬ 
namics,  Without  the  gimbals,  all  of  the  wheels  always  begin  unstable 
whirling  inmediately  above  the  first  critical  speed  (typically  about  50- 
100  ri^m) .  In  addition,  it  was  found  that  whirl  stability  is  sensitive  to 
the  support  rope  configuration.  Die  most  stable  configuration  was  found  to 
be  tivisted  ("maypole")  support  ropes,  approximating  a  synthetic  quill  shaft. 

Die  most  impressive  feature  of  the  flexible  fly\>:heel  demonstrated  by 
the  experiments  was  self-balancing,  a  result  of  the  flexible  supports. 

None  of  the  fl)nvheels  had  to  be  balanced,  and  no  noticeable  synchronous 
vibration  was  ever  transmitted  to  the  structure. 

Maximum  speeds  achieved  with  the  three  wheels  were  2600  rpm  (10-1/2" 
wheel),  2000  rpm  (22-1/2"  wheel)  and  1500  rpm  (30"  wheel).  With  the  twisted 
support  ropes,  these  speed  limits  were  imposed  by  electrical  power  supply 
limitations  in  overcoming  air  drag.  With  the  untwisted  ropes,  speeds  were 
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limited  by  the  subsynchronous  whirl  instability  threshold. 

A  computerized  stability  analysis  was  also  developed  to  guide  the  design., 
improve  understanding  o£  the  destabilizing  mechanism,  and  ultimately  to 
optimize  the  design  parameters  for  dynamic  stability.  All  of  the  predicted 
critical  speeds  and  mode  shapes  agreed  closely  with  experimental  measurements , 
and  the  thresholds  of  stability  could  be  accurately  simulated  by  adjusting 
the  (unknown)  values  of  internal  hysteresis  in  the  mathematical  model.  One 
result  of  the  analysis  is  a  theory  to  rationally  explain  why  the  maypole  con¬ 
figuration  is  stable  and  the  untwisted  support  rope  configuration  is  not. 

In  order  to  answer  questions  about  the  torque  capacity  of  the  twisted 
maypole  support  rope ,  measurements  of  the  actual  air  drag  torque  on  the 
22  1/2"  flywheel  were  made  at  various  speeds  which  showed  a  power  capability 
to  generate  more  than  1.3  lOV  at  6,000  rpm,  and  more  than  2.2  KW  at  10,000 
rpm.  On  several  occasions  during  the  experiments,  the  electrical  current 
demand  due  to  air  drag  caused  circuit  breakers  to  trip  in  the  power  supply, 
thus  subjecting  the  flywheel  to  a  severe  transient  decelerating  torque.  It 
was  found  that  the  flywheel  could  handle  these  transients  with  no 
unfavorable  effects. 

The  production  cost  for  a  conceptual  50  KWh  residential  flexible 
flywheel  energy  storage  system,  is  estimated  to  be  $24,160  ($483/KlAJh) .  Based 
on  that  cost  estimate  and  using  a  life  cycle  cost  methodology  developed  for 
solar  energy  systems,  the  residential  present  value  cost  ($1979)  for  this 
50  KWh  system  would  be  $36,437. 

Probably  the  most  important  conclusions  to  be  drawn  from  this  research 

are : 
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(1)  Subsynchronous  whirl  instability  at  supercritical  speeds  due  to 
internal  friction  is  strongly  suppressed  by  a  giniballed  motor/ 
generator  support  system.  (ITiis  finding  is  applicable  to  more 
conventional  flyivheels,  as  well  as  to  the  flexible  flywheel.) 

(2)  A  flexible  flyivheel  made  from  synthetic  rope  or  fiber  with  no 
bonding  agent  is  easy  to  construct  and  can  be  made  self¬ 
balancing  by  a  properly  designed  support  system. 

(3)  The  flexible  flywheel  system  offers  a  relatively  inexpensive 
and  efficient  tool  for  investigating  flywheel  rotor  dynamics  in 
general,  since: 

(a)  failures  are  not  catastrophic, 

(b)  critical  speeds  are  so  low  that  whirl  modes  can  be 
visually  observed,  and 

(c)  support  rope  configurations  can  be  varied  to  approxi¬ 
mate  the  shaft  constraints  of  conventional  flywheels. 

(4)  The  simplicity  of  the  flexible  flywheel  construction  promotes 
economical  production  using  conventional  manufacturing 
processes  and  equipment. 


CONFIGURATION  AND  ADVANTAGES  OF  THE  FLEXIBLE  FLYWHEEL 


Figure  1  shows  a  photograph  o£  an  early  experimental  flexible  flywheel. 
Notice  that  the  support  ropes  carry  only  the  weight  of  the  flywheel  and 
that  there  is  no  rigid  hub  or  spoke  array  as  in  a  conventional  solid  or  com¬ 
posite  flywheel. 

The  radial  ropes  constraining  the  flywheel  to  the  lower  circular  flange 
were  intended  to  control  whirling,  an  approach  which  was  subsequently  found 
to  aggravate  the  problem.  In  the  present  design,  radial  stiffness  of  the 
rotor  is  provided  by  gravity  only,  as  in  a  pendulum. 

The  flywheel  is  a  simple  hoop,  constructed  by  coiling  ropes  or  fibers 
in  a  circle  of  the  desired  radius  until  the  design  mass  is  achieved.  At 
speed,  centrifugal  force  maintains  the  circular  configuration.  For  the  sake 
of  shape  integrity  when  not  spinning,  one  layer  of  rope  is  wound  around  the 
cross  section  of  the  hoop.  (Early  models  were  of  macrame  construction,  but 
this  was  found  to  sacrifice  strength  with  no  apparent  advantage.)  The 
loose  ends  of  the  coiled  hoop  are  joined  by  tying  or  splicing. 

From  the  standpoint  of  energy  density,  this  fiber/hoop  configuration  is 
the  optimum  attainable.  It  is  obviously  cheap  to  manufacture,  and  newly 
developed  high  strength  fibers  will  be  easy  to  incorporate  into  the  design 
without  the  need  to  develop  new  bonding  techniques. 

The  question  of  torque  capacity  and  acceleration  has  been  investigated 
for  the  support  rope  configuration  shown  in  Figure  2  and  appears  not  to  be 
a  problem.  The  most  stable  configuration  (from  the  standpoint  of  whirling) 
occurs  when  the  support  ropes  are  twisted  around  each  other.  The  torque 
capacity  of  this  configuration  has  also  been  investigated  and  is  ample  for 
the  high  design  speeds  achievable  with  synthetic  fibers. 
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Flexible  Flywheel  Without  Lower  Radial-Constraint  Ropes 
(Flywheel  Shown  in  Sub  synchronous  t^^hirl  Mode) 


The  advantages  claimed  for  the  flexible  flywheel  with  proper  support 
system  are: 

1.  High  strength  fibers  in  pure  tension  with  no  bonding  materials  to 
create  mismatch  in  elasticity  or  strength. 

2.  Self-balancing,  due  to  a  highly  flexible  rotor  operating  at  super¬ 
critical  speeds , 

3.  Simple  construction  promotes  economical  production.  Newly 
developed  fibers  can  be  easily  incorporated  without  expensive 
redevelopment  of  manufacturing  techniques. 

4.  Unbonded  fiber  construction  simplifies  the  safe  containment  of  a 
rotor  failure. 

5.  Large  hoop  diameters  made  possible  by  the  simple  construction 
allow  low  operating  speeds  and  attendant  reduced  power  losses. 

FLEXIBLE  FLYtVHEEL  DESIGN  EQUATIONS  AND  CONSTRAINTS 

Analysis  has  shown  that  to  optimize  energy  density,  a  flywheel  should 
be  constructed  of  high  strength  materials  to  operate  at  high  speeds .  The 
maximum  storage  energy  in  a  hoop  flywheel  is 

E  =  TrrP  ft-lb 

E=  (.3768)  (10"^)7rRP  KW-hrs  (1) 

where  R  -  hoop  mean  radius,  ft 

P  -  cumulative  strength  of  all  hoop  fibers,  lb. 

Notice  that  the  material  mass  density  does  not  appear  in  the  equation. 
Materials  with  high  mass  density  (heavy  materials)  do  not  optimize  energy 
density,  which  is  contrar>"  to  popular  intuition. 

The  mass  density  does,  however,  affect  the  speed  at  which  the  maicimum 

energy  is  stored.  This  speed  is  given  by 


18 


Vvhere  w  =  hoop  specific  weight,  lb/£t 
g  =  32.2  ft/sec2 

Although  there  aie  many  advantages  to  using  super-strong  fibers  to  take 

advantage  of  equation  (1) ,  the  resulting  high  rotational  speeds  (most  high- 

strength  fibers  are  not  heavy)  pose  rotor  dynamics  and  bearings  problems 

which  must  be  properly  appreciated  in  the  preliminary  design  phases  of  any 

modem  flywheel.  For  example,  contemporary  electrical  motors  and  generators 

are  designed  to  operate  well  below  the  speed  dictated  by  equation  (2)  for  a 
R 

Kevlar  wheel. 

There  are  also  practical  constraints  on  the  dimensions  of  the  flywheel  to 
be  used  in  the  home  or  on  a  small  farm.  Figure  3  gives  dimension  limits  for  a 
fiber  hoop  flywheel. 

Applying  the  above  equations  and  constraints  to  the  design  of  a  flexible 

flywheel  allows  the  calculation  of  hoop  sizes,  weights  and  costs  for  various 

R  R 

choices  of  material.  Comparative  results  for  Dacron  ,  steel,  and  Kevlar  are 

shown  in  Table  1,  for  a  6  ft.  diameter  10  KWh  flywheel.  Values  for  a  6  ft. 

diameter  50  KWh  wheel  are  given  in  Table  4,  page  59. 

For  a  successful  energy  storage  system,  reasonable  cost  is  essential. 

R 

Table  1  shows  that  the  high  strength  of  Kevlar  pays  off  in  a  lower  total 
material  cost  for  the  flywheel.  The  technical  price,  however,  is  high  oper¬ 
ating  speeds.  It  will  be  shown  in  the  next  section  that  the  speed  can  also 
be  brought  down  by  the  use  of  a  larger  hoop  diameter.  There  are  many  appli¬ 
cations,  however,  where  space  is  limited  and  Figure  3  or  similar  constraints 
will  govern  the  design. 
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Table  1.  Parametric  Values  for  a  6  ft.  Diameter  10  KWh  Flexible  Flywheel 

(Safety  Factor  =2) 


MATERIAL 

HOOP  l^IGHT 
lbs. 

MAX.  RPM 

LENGTH  L 
ft. 

MAT'L  COST 
$ 

Dacron 

(1/4") 

1298 

5,167 

1.74 

2954 

Steel 

(1/2"IRWC) 

2123 

4,039 

.642 

1938 

Kevlar 

(1500  Den. -"29") 

162 

14,606 

.287 

1280 

r  ^  4  too 


L  ^  4  fto 


Fig.  3  Dimension  Constraints 
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THE  EFFECT  OF  SIZE 

There  are  significant  advantages  to  increasing  the  hoop  diameter  of  a 
flexible  flyTAfheel. 

The  values  in  Table  1  were  calculated  under  the  constraints  of  Figure 
3,  If  we  choose  Kevlar^  as  the  material,  remove  the  constraint  on  the  hoop 
radius  ”R”,  and  again  use  a  safety  factor  of  2,  equations  (1)  and  (2)  can  be 
written  as 


(10.97)R  KWh 

CD’ 

and 

=  31,005/R  rpm 

(2)’ 

lliese  equations  show  that  tlie  energy  stored  is  increased  and  operating 
speed  is  decreased  by  making  the  hoop  diameter  larger  (with  constant 
circular  cross  section) .  Reference  4  shows  that  the  latter  effect  (lower 
speed)  reduces  power  losses. 

D 

Thus,  a  10'  diameter  Kevlar  flexible  flywheel  could  store  50  KWh  at 
a  speed  of  only  6200  rpm,  and  have  bearing  power  losses  of  only  178  watts. 

(See  Appendix  A  for  power  loss  calculation  method.) 

Furtliermore ,  a  flexible  flywheel  of  large  diameter  can  be  constructed  with¬ 
out  the  problems  which  would  be  expected  from  composite  construction  (inhomo¬ 
geneous  curing,  uncertain  shape  and  balance  retention,  high  cost  of  manufacture, 
etc.) . 

ROTOR  DYNAMICS 

It  was  recognized  early  that  the  significant  technical  problems  associated 
with  the  flexible  flywheel  would  be  in  the  area  of  rotor  dynamics,  specifically 
the  problem  of  subsynchronous  whirling  due  to  internal  friction.  The  advantage 
of  self-balanced  operation  at  high  supercritical  speeds  must  be  purchased  with 
the  price  of  suppressing  or  ax'^oiding  a  self-excited  dynamic  instability.  This 
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is  a  challenge  with  a  high  payoff  for  success,  and  one  which  the  author  believes 
has  been  met. 

Ever  since  Jeffcott's  analysis  of  synchronous  rotor  whirl  in  1919,  rotor 
dynamicists  have  known  that  a  flexible  rotor  displays  a  "critical  speed  inver¬ 
sion",  in  which  the  center  of  mass  comes  inside  the  whirling  rotor  center  at 
sui^ercritical  speeds. 

Tlie  flexible  flywheel,  by  virtue  of  its  low  stiffness  rotor,  always  operates 
at  speeds  which  are  highly  supercritical,  where  Jeffcott  showed  that  synchronous 
whirl  (due  to  unbalance)  amplitudes  are  minimized.  Experiments  to  date  have 
verified  that  the  flexible  flyivheel  produces  extremely  low  levels  of  synchronous 
vibration,  with  no  precision  balancing  required. 

Not  long  after  Jeffcott's  results  became  widely  known  and  applied,  it  was 
learned  that  rotating  machinery  can  become  dynamically  unstable  in  subsynchronous 
\shirl  at  supercritical  speeds,  if  the  ratio  of  internal  friction  (in  rotating 
parts)  to  external  damping  is  high  enough. 

Early  tests  of  the  flexible  flyv/heel  showed  subsynchronous  whirling  iv'-hich 
tended  to  grow  with  speed  and/or  time.  Figure  4  illustrates  the  mechanism  of 
the  internal  friction  excitation.  For  subsynchronous  whirl,  the  spin  speed  0. 
is  faster  than  the  whirl  speed  tj).  As  support  rope  3  moves  around  to  position 

1,  its  rate  of  strain  is  a  maximum  at  position  2,  thus  generating  the  friction 

force  F  or  the  hoop  which  is  tangential  to  the  whirl  orbit  in  the  forward  direction. 

Rotor  dynamics  theor}''  and  analysis  has  identified  several  ways  of  suppressing 
self-excited  subsynchronous  whirl.  They  are; 

1.  Flexible  bearing  supports. 

2.  .4s)nnmetric  bearing  support  stiffness. 

3.  Bearing  support  damping. 

4.  Bearing  support  mass  (d>Tiamic  absorber  effect). 
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Fig.  4.  Instantaneous  Configuration  of  Rope  Ring 

Before  describing  how  these  methods  were  applied  to  stabilizing  a  flexible 
flwheel,  it  is  useful  to  take  a  preliminary  look  at  the  other  "system”  design 
considerations.  A  more  detailed  system  design  and  cost  analysis  is  presented 
in  the  last  section. 

PRELIMINARY  SYSTEM  DESIGN  CONSIDERATIONS 

1.  A  shaft  seal  through  the  vacuum  chamber  wall  is  expensive. 

2.  A  disconnect  clutch  also  increases  the  total  system  cost. 

3.  A  new  motor/generator  must  be  developed  to  match  fly\\^heel  torque- 
speed  characteristics. 

4.  Low  friction  bearings  must  be  designed  for  the  application,  to 
operate  in  a  vacuum  environment. 

5.  The  total  number  of  bearings  should  be  minimized,  for  lowest  cost, 
and  for  lowest  power  loss. 

IVhen  these  design  considerations  are  coupled  with  the  design  requirements 
to  suppress  the  subsynchronous  whirl,  a  design  philosophy  for  the  flexible  fly¬ 
wheel  emerges.  Tlie  basic  elements  of  this  philosophy  are  shown  in  Table  2. 
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Table  2.  Flexible  Flywheel  Design  Philosophy 


Design  Factors  or  Constraints 
Shaft  seal  is  expensive 
New  motor/generator  required 


Solution  or  Conclus ion 


Put  motor  inside  vacuum  chamber 


Clutch  increases  cost 
New  bearings  required 
fhnimize  no.  of  bearings 


Support  flywheel  directly  from 
motor  shaft 


Need  low  support  stiffness 
Need  asymmetric  support  stiffness 
Need  support  damping  and  mass 


Gimbal  motor/ generator  on 
nonintersecting  axes. 


THE  GIMBALLED  SUPPORT  DESIGN  MD  E^CPERIMENTAL  RESULTS  OF  IVHIRL  STABILITY 

Figure  5  shows  how  the  motor/generator  can  be  gimballed  on  nonintersecting 
axes  to  provide  the  low  support  stiffness,  stiffness  asymmetry,  and  bearing 
support  mass  (the  motor  itself) ,  which  are  the  parameters  inportant  to  whirl 
stability. 

The  motor  bearings  are  designed  to  support  the  flywheel,  thus  minimizing 
the  number  of  bearings  and  eliminating  the  necessity  for  a  clutch,  since  the 
motor  aimature  becomes  effectively  a  part  of  the  fl)^Aeel  inertia. 

Work  on  this  contract  has  been  aimed  at  verifying  whirl  stability  for  this 
concept,  both  ex])eriraentally  and  analytically.  Tliis  section  describes  the  ex¬ 
perimental  work;  analysis  is  described  in  the  section  following. 
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Fig.  5  Flexible  Flyulieel  On  Gimballed  Support 


A  small  scale  model  was  constructed  for  preliminary  evaluation  of  whirl 
stability  characteristics.  The  distance  h  from  the  Z  axis  down  to  the  motor 
center  of  mass  was  made  adjustable,  so  as  to  vary  the  stiffness  asymmetry 
produced  by  the  gimbals.  The  flywheel  used  in  this  model  is  10  1/2”  diameter 


and  is  made  from  a  continuous  coil  of  Nylon  rope.  Th,e  weight  of  the  fl>n^heel 
hoop  alone  is  2.7  lb.  Figures  6  and  7  are  photographs  of  the  small  model. 

Tests  were  made  with  the  gimbals  locked  (rigid) ,  and  with  the  gimbals 
free,  to  evaluate  the  effect  on  whirl  stability. 

Nitli  locked  gimbals,  the  small  flywheel  is  violently  unstable  in  subsyndir 
nous  whirl  at  all  speeds  above  100  rpm,  for  all  support  rope  configurations. 

With  the  gimbals  free,  and  twisted  "maypole"  support  ropes,  the  flywheel 
is  stable  and  runs  smoothly  up  to  the  maximum  speed  allowed  by  available  motor 
power  to  overcome  air  drag  (about  2600  rpm) .  With  the  gimbals  free  and  un¬ 
twisted  support  ropes,  the  threshold  speed  of  instability  is  between  600-750 
rpm,  depending  on  the  length  of  the  support  ropes.  In  this  configuration, 
asymmetry  of  the  gimbal  supports  (one  gimbal  axis  free,  one  locked)  produces 
a  slightly  higher  threshold  speed  than  with  both  gimbals  free. 

No  effort  has  ever  been  made  to  balance  this  flywheel,  and  at  one  time 
small  unbalance  weights  were  added  at  random  to  the  twisted  rope  configuration 
to  demonstrate  its  insensitivity  to  unbalance.  The  smooth  and  stable  operation 
was  not  impaired. 

A  larger  gimbal- supported  motor/ generator  system  was  used  to  test  a  22-1/2 
diameter  Dacron  flywheel  weighing  19  lb.  Its  15  HP  motor  allowed  speeds  in  the 
atmosphere  up  to  2000  rpm  before  air  drag  caused  power  supply  circuit  breakers 
to  trip. 

Figures  8,  9,  and  10  show  the  22-1/2"  wheel  with  its  motor/generator, 

gimbal  support  system,  and  protective  cage. 

The  major  shortcoming  of  the  22-1/2"  wheel  was  its  relatively  low  weight 
(19  lb.)  compared  to  the  motor/ generator  weight  (90  lb.).  It  was  felt  that 
this  impeded  the  stabilizing  influence  of  the  gimbals,  an  intuition  which  was 
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Gimballed  Motor  For  Small  Model  (10-1/2")  Flywheel,  Top 
and  Side  Views 
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Fig.  9.  22-112"  Flyivheel 


later  verified  bv  results  frorn  the  coinputer  stability  analysis.  Experiments 
with  the  22-112"  wheel  showed  an  even  stronger  sensitivity  to  the  method  of 
support  rope  attachment,  both  at  the  flywheel  and  at  the  motor  hub.  The  2,000 
rpm  speed  was  achieved  with  untwisted  steel  cable  supports  attached  to  the  fly¬ 
wheel  with  rope  bridles.  As  mentioned  above,  the  speed  limitation  in  this  case 
was  imposed  by  excessive  current  demand  tripping  circuit  breakers  in  the  power 
supply.  Shortly  after  this  test,  the  22-1/2”  wheel  was  removed  so  that  the 
motor/generator  and  gimbal  supports  could  be  installed  in  a  vacuum  chamber 
(as  a  safety  measure,  and  for  possible  pump-down  at  a  later  date).  After  this 
was  done,  the  22-1/2”  wheel  was  reinstalled,  but  the  2,000  rpm  speed  could  not 
be  repeated  with  untxvisted  support  ropes.  Since  a  considerable  volume  of  air 
blows  out  of  one  of  the  chamber  ports  when  the  flyiAeel  runs,  it  is  felt  that 
aerodynamic  pressure  may  produce  an  additional  destabilizing  force  in  the  chamber. 
Another  possibility  is  that  the  support  cable  attachment  scheme  was  not  precisely 
duplicated  after  re installation. 

Finally,  a  30”  diameter  Nylon  wheel  weighing  88  lb.  was  constructed  and 
tested  in  the  vacuum  chamber,  using  the  15  HP  motor/gimbal  assembly  shown  in 
Figure  8.  This  flywheel  can  store  1  KlV-hr  at  10,650  rpm,  and  has  a  breaking 
speed  of  14,548  rpm.  Figure  11  shows  the  vacuum  chamber,  and  Figure  12  is  a 
photograph  of  the  30”  flyi\'heel  in  the  vacuum  chamber  (looking  dovm  from  the  top). 
With  the  flywheel  supported  by  untwisted  support  ropes  made  from  monofilament 
Nylon,  the  threshold  speed  of  instability  was  1,000  rpm.  To  simulate  the  twisted 
”ma>^ole”  support  ropes  which  had  proved  stable  in  the  smaller  models,  and  yet 
retain  torque  capability  in  both  directions,  a  four  strand  bridle  was  braided 
into  a  square  cross  section.  With  this  support  rope  configuration,  the  maximum 
speed  achieved  was  1500  rpm  and  was  limited  only  by  the  power  supply  overload 
due  to  air  drag. 
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TORQUE  AND  POWER  CAPABILITY  OF  THE  SUPPORT  ROPES 


To  answer  questions  raised  about  the  ability  of  the  support  ropes  to 
transmit  sufficient  torque  for  generating  useful  magnitudes  of  electrical 
power,  the  torque  actually  transmitted  in  overcoming  air  drag  on  the  22  1/2” 
flywheel  was  measured.  This  was  accomplished  by  measuring  the  rope  angles 
in  photographs  and  movies,  and  relating  these  angles  to  static  torque 
measurements.  The  results,  converted  to  KW  at  various  speeds,  are  shown  on 
Figure  13.  It  can  be  seen  that  the  twisted  rope  configuration  (previous 
Figure  10),  the  one  of  primary  interest,  would  generate  over  1.3  KIV  at 
6,000  rpm  and  over  2.2  KIV  at  10,000  rpm.  It  can  be  assumed  that  this 
torque/power  transmitting  capability  will  increase  with  the  flywheel  diame¬ 
ter  and  with  the  moment  of  inertia  of  the  cross-section  area  of  the  quill 
shaft  which  is  created  by  the  twisted  support  ropes.  This  assumption  is 
supported  by  torque  tests  on  the  30  inch  diameter  nylon  flywheel  in  which 
four  1/4  inch  support  ropes  (same  size  ropes  as  used  on  the  22  1/2  inch 
diameter  flywheel)  were  square-braided  to  form  a  quill  shaft  with  a  cross- 
section  area  of  greater  moment  of  inertia  than  was  obtainable  by 
"maypoling”  the  ropes  into  a  shaft  of  essentially  circular  cross-section. 

The  input  power  required  for  the  motor  to  drive  this  30  inch  diameter  fly¬ 
wheel  at  1500  rpm  in  the  atmosphere  translates  into  approximately  6  KIV 
power  generating  capability  in  a  vacuum, 

INSTRUMF:ITrATION  AND  MEASUREMEOTS 

Flywheel  speed  of  the  small  (10-1/2”)  model  was  measured  with  a  photocell 
reflective  tachometer.  The  reflective  tape  was  installed  on  the  flywheel  hoop 
itself.  Natural  frequencies  were  measured  with  a  SELSPOT  optical  tracking  in¬ 
strument,  tracking  the  motion  of  LED's  installed  on  the  flywheel.  The  signal 
xvas  captured  on  a  storage  oscilloscope,  so  that  the  logarithmic  decrement 
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Fig.  13.  Electrical  Power  Capability  of  Flexible  Flywheel  From 
Measured  Torque 

fdiimping)  could  be  determined  as  tvell  as  the  natural  £requenc)w 

Flwheel  speeds  of  the  larger  (22-1/2”  and  30”)  flywheels  were  measured 
with  a  pulse  tachometer,  consisting  of  a  magnetic  transducer  excited  once  per 
revolution  by  a  projection  from  the  steel  hub  of  tlie  motor  shart,  with  the 
period  between  pulses  measured  by  an  H-P  electronic  timer.  (M  electronic 
tachometer  reading  in  ipm  was  also  used  in  the  early  stages  of  the  project) . 

UTiirling  frequencies  of  instability  were  so  low  (typically  1-2  Hz)  that 
it  was  found  possible  to  measure  them  visually,  using  a  stopwatch,  and  by 
movies,  placing  a  clock  with  a  sweep  second  hand  in  the  camera  field  of  view'. 

Due  to  the  flexible  sipport  with  no  hub  or  spokes,  the  £ly\fheel  hoop 
whirled  as  a  rigid  body.  No  flexural  vibrations  of  the  hoop  (as  envisioned 
in  the  original  proposal  for  this  contract)  were  ever  observed. 


36 


FLEXIBLE  FLYIVHEEL  STABILITY  A^IALYSIS 


I.  The  Mathematical  Model 

A  mathematical  stability  analysis  was  performed  on  the  flexible  fly¬ 
wheel  so  that  its  stability  characteristics  can  be  predicted  for  any  given 
set  of  system  parameters.  In  particular,  it  is  desirable  to  know  which 
parameters  have  the  most  pronounced  effect  on  the  system  stability  when  they 
are  changed.  Three  different  mathematical  models  were  considered  and  actual 
computer  simulation  was  performed  for  two  of  them.  The  most  important  dif¬ 
ference  between  the  three  models  is  the  configuration  of  the  support  ropes 
which  join  the  flywheel  to  the  gimballed  motor. 

A.  Rigid  shaft  model 

The  first  model  to  be  considered  is  sho™  in  Figure  14.  The  gimballed 
support  system  has  been  modeled  to  allow  for  the  rotation  of  the  motor  about 
the  gimbal  axes  in  the  presence  of  viscous  damping.  All  mass  properties  of 
the  motor  and  its  support  system  have  been  taken  into  account ,  including  the 
gyroscopic  effects  of  the  motor's  rotating  parts.  The  flywheel  itself  is 
modeled  as  a  rigid  rotor  with  all  its  corresponding  mass  properties. including 
gyroscopics.  External  damping  on  the  flywheel  is  assumed  to  be  viscous  air 
drag  proportional  to  the  flywheel's  translational  velocity  and  acts  at  the 
flywheel's  center  of  mass. 

The  most  important  aspect  of  this  model  is  the  modeling  of  the  support 
ropes.  The  support  ropes  are  assumed  to  be  straight  and  to  remain  straight 
at  all  operating  speeds.  This  would  result  in  the  ropes  behaving  collective¬ 
ly  as  a  rigid  shaft  attached  rigidly  to  the  flywheel  and  flexibly  to  the 
motor.  With  the  ropes  modeled  as  such,  the  internal  friction  force  due  to 
hysteresis  in  the  ropes  will  then  be  proportional  to  the  misalignment  angle 
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i 

f 

i  between  the  motor *s  main  axis  and  the  rigid  shaft  that  makes  up  the  support 

t 

j  ropes . 

^  The  model  just  described  has  five  degrees  of  freedom.  Two  angles  a  and 

'  3  have  been  used  to  describe  the  motor's  orientation,  each  angle  being  the 

j  rotation  about  one  of  the  gimbal  axes.  Three  more  angles,  (j> ,  0  and  are  used 

I 

!  to  specify  the  flywheel's  orientation.  One  of  these  flywheel  angles,  i};,  con- 

I 

I  stitutes  the  spin  of  the  fly\dieel,  and  by  assuming  the  flyw^heel's  acceleration 

( 

I  to  be  small,  constant  speed  can  be  assumed  to  reduce  the  system  to  four  degrees 

i 

<  of  freedom. 

;  Since  the  model  has  four  degrees  of  freedom,  the  system  will  have  four 

i 

!  natural  frequencies  and  four  corresponding  mode  shapes .  Each  mode  shape  is , 

!  in  fact,  a  whirling  configuration  with  both  the  flywheel  and  the  motor  executing 

/ 

■  a  circular  or  elliptical  path.  Two  of  these  modes  will  have  the  flywheel  exe- 

t 

! 

;  cuting  a  rather  large  motion  with  the  motor  following  along  with  it.  One  of 

i  these  will  be  a  forward  whirl  in  the  direction  of  flywheel  rotation  and  one 

•  will  be  backwards  whirl.  The  two  remaining  mode  shapes  occur  at  the  critical 

f 

i  speeds  of  the  motor  support  system.  In  these  modes  the  motor  will  execute  the 

i  large  motion  and  the  flywheel  will  remain  relatively  stationary.  Once  again, 

'  one  mode  will  be  forward  whirl  and  the  other  will  be  backwards. 

The  four  second  order  differential  equations  of  motion  for  this  system 
have  been  derived  using  Lagrange's  method  and  are  given  in  Appendix  B. 

B.  Bowed- out  rope  model 

blien  the  flywheel  is  spinning  at  a  high  speed,  the  support  ropes  no  longer 
remain  straight  but  will  bow  outwards  due  to  the  centrifugal  forces  acting  on 
them,  see  Figure  15.  Witii  the  support  ropes  in  this  bowed-out  configuration, 
they  no  longer  act  like  a  rigid  shaft,  rigidly  connected  to  the  flywheel  as  was 
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assumed  in  the  first  model.  Tliis  additional  flexibility  provided  by  the 
bowed-out  ropes  allows  the  flywheel  to  take  on  any  orientation  independent 
of  its  position. 

Mien  the  flywheel  is  spinning  at  a  high  speed,  the  gyroscopic  moments 
tending  to  keep  the  flywheel  horizontal  are  very  strong.  This  gyroscopic 
effect  tends  to  keep  the  flywheel  horizontal,  even  when  it  is  in  a  whirling 
configuration.  Tliis  causes  the  internal  friction  forces  to  act  differently 
than  in  the  rigid  shaft  mode.  The  internal  friction  force  in  each  rope  is 
proportional  to  the  rate  of  strain  which  in  turn  is  determined  by  the  geometry 
of  the  system  and  the  speed  of  the  flywheel.  But,  since  the  geometr)'’  is  now 
different,  then  so  will  be  the  internal  friction  force  at  any  given  speed. 

Taking  into  account  the  new  geometry  of  the  system,  new  differential 
equations  of  motion  were  derived  again  using  Lagrange's  method,  see  Appendix  B, 
Tlie  same  angular  coordinates  were  used  as  for  the  rigid  shaft  model,  and  are 
analogous  to  their  counterparts  for  the  rigid  shaft  model. 

C.  Twisted  rope  model 

The  third  model  considered  is  essentially  the  same  as  the  first  two  models. 
The  main  difference  is  that  the  support  ropes  are  now  twisted  around  them¬ 
selves,  see  Figure  16.  This  twisting  of  the  ropes  causes  a  considerable  change 
in  the  stability  of  the  system,  as  is  observed  with  the  small  (10-1/2”)  model 
flywheel. 

For  this  twisted  rope  system  the  differential  equations  of  motion  were  not 
derived  and  a  computer  analysis  was  not  performed.  However,  a  theory  has  been 
developed  to  explain  the  observed  behavior  of  the  twisted  rope  system.  (Both 
its  stability  with  the  gimballed  motor  and  its  instability  with  the  rigidly 
supported  motor.) 
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Fig.  16.  Flexible  Flyu'heel  ?lodel  on  Twisted  Support  Ropes 
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Mien  the  support  ropes  are  tivisted  and  the  flywheel  is  in  a  whirling 
configuration,  see  fig  16,  gyroscopic  moments  will  keep  the  flywheel  hori¬ 
zontal  as  with  the  bowed-out  rope  model .  This  produces  internal  friction 
forces  in  four  different  regions  of  the  support  ropes,  labeled  A,  B,  C  and 
D  in  the  figure.  In  regions  A  and  D,  the  internal  friction  forces  are  pro¬ 
duced  in  exactly  the  same  manner  as  for  the  rigid  shaft  model  and  thus 
will  be  in  the  direction  of  forward  whirl.  In  region  B  the  internal  fric¬ 
tion  force  is  produced  in  mudi  the  same  manner  as  in  a  spline  fit  coupling. 

It  will  therefore  be  in  the  direction  of  forward  whirl  and  it  will  be  equal 
to  zero  when  the  t^^risted  ropes  are  in  line  with  the  motor.  The  internal 
friction  forces  produced  in  region  C  are  produced  in  the  same  way  as  for 
region  B  except  that  now  it  will  be  in  the  direction  of  backwards  whirl. 

Mnen  the  rotating  flywheel  is  in  a  whirling  configuration  with  the 
motor  held  rigid,  tiae  forces  produced  in  regions  B  and  C  will  be  equal  and 
opposite  and  thus  will  cancel.  The  forces  produced  in  regions  A  and  D  avill 
both  be  trying  to  produce  fonvard  whirl ,  and  tlius  the  flywheel  should  be 
expected  to  go  unstable  at  a  rather  low  speed.  This  prediction  matches 
exactly  the  behavior  observed  for  the  small  (10-1/2”)  model  flyivheel. 

hlien  the  motor  is  allowed  to  rotate  in  its  gimballed  support,  and  the 
flwheel  is  in  a  whirling  configuration,  the  motor  will  follow  the  flywheel 
around  so  that  the  motor  axis  is  always  in  line  with  the  twisted  ropes. 

This  action  will  cause  the  friction  forces  in  regions  A  and  B  to  both  be 
zero.  Tlie  friction  force  in  region  C  will  still  be  in  the  direction  of 
backwards  whirl,  and  in  region  D  it  will  be  forward.  Thus,  the  only  internal 
friction  forces  present  will  be  acting  against  each  other.  There  is  no  rea¬ 
son  to  suspect  that  they  will  be  equal  and  opposite,  but  the  resultant  can 
be  expected  to  be  much  smaller  than  in  the  case  of  the  rigid  shaft  model  or 
the  bowed-out  ropes  model.  It  also  should  he  possible  to  somehow  adjust  the 
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two  forces  so  that  they  are  equal  and  opposite,  resulting  in  no  net  internal 
friction  force  and  stable  operation  at  all  speeds. 

Ihe  small  model  flywheel  is  seen  to  be  stable  at  all  speeds  of  operation 
when  the  support  ropes  are  twisted  and  the  motor  free,  and  so  the  predictions 
based  on  the  theory  just  presented  match  exactly  its  observed  behavior.  If 
this  theory  should  prove  to  be  correct,  then  the  internal  friction  instability 
can  be  eliminated  even  if  the  flywheel  is  operated  in  a  vacuum  without  the 
stabilizing  effect  of  air  drag. 

II.  Internal  Friction  Force 

In  both  the  rigid  shaft  model  and  the  bowed-out  ropes  model,  the  internal 
friction  force  is  generated  in  the  same  manner  by  hysteresis  in  the  support 
ropes.  The  magnitude  and  direction  of  the  internal  friction  force  is  dependent 
upon  the  displacement  of  the  system  although  this  dependence  is  entirely  differ¬ 
ent  for  the  two  models.  Hie  qualitative  relationship  between  the  system  dis¬ 
placement  and  the  direction  of  the  internal  friction  force  will  now  be  considered 
for  the  two  models. 

A.  Rigid  shaft  model 

For  this  model,  the  internal  friction  force  is  proportional  to  any  mis¬ 
alignment  between  the  motor  and  the  flyvdieel.  The  internal  friction  force  will 
be  zero  only  if  the  axis  of  rotation  of  the  flywheel  is  colinear  with  the  axis 
of  rotation  of  the  motor.  Figures  17a  and  17b  show,  in  a  simple  way,  all  the 
possible  configurations  that  the  system  can  assume.  The  motor  is  shown  in  an 
arbitrarily  deflected  position  and  the  flywheel  can  be  anywhere  on  the  indicated 
curve  for  the  given  motor  deflection.  Depending  on  where  the  flwheel  is  on  this 
curve,  the  friction  force  may  be  in  the  direction  of  forv'^ard  whirl  or  backward 
whirl.  The  various  regions  are  also  indicated  on  the  figure. 
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Fig.  17a.  Internal  Friction  Force  For  Model  With  "Rigid- Shaft” 
Support  Ropes  For  Deflection  About  Upper  Gimbal  Axis 


set  of  all 
possible  positions 
for  flyv.'lieel  center 


position  of  zero 
friction  force 


'b.  Internal  Friction  Force  For  Model  With  "Rigid  Shaft" 
Support  Ropes  For  Deflection  About  Lower  Gimbal  Axis 


If  a  mode  shape  which  executes  foru^ard  whirl  should  have  the  flywlieel 
in  a  region  labeled  "forward”  on  the  figure,  then  the  internal  friction 
force  will  be  in  the  direction  that  the  flywheel  is  trying  to  move.  Under 
these  circumstances,  the  system  will  stand  a  very  good  chance  of  going  un¬ 
stable,  Tlie  same  holds  true  if  a  backwards  mode  shape  should  happen  to  put 
the  flywheel  in  a  region  labeled  "backwards". 

On  the  other  hand,  should  a  forv\'ard  mode  shape  put  the  flywheel  in  a 
region  labeled  "backvAfards",  tlien  the  flywheel  would  not  be  expected  to  go 
unstable  since  the  internal  friction  force  would  be  opposite  the  direction 
that  the  fl>avlieel  wants  to  go.  A  backivards  mode  shape  in  a  "forvard"  region 
would  produce  the  same  effect. 

Conceivably,  it  could  be  possible  to  find  a  special  set  of  system  para¬ 
meters  that  will,  for  each  mode  shape,  locate  the  flyv\rheel  in  a  region  of 
opposing  internal  friction  force.  If  this  is  done,  then  one  can  expect  stable 
operation  at  any  operating  speed  for  a  flywheel  that  actually  is  supported  by 
a  rigid  shaft,  or  quill  shaft. 

B.  Bowed- out  rope  model 

For  this  model,  the  relationship  between  system  deflection  and  the  in¬ 
ternal  friction  force  is  more  complicated  and  depends  on  more  system  para¬ 
meters.  Uie  tvTO  drawings  in  Figure  IS  show  the  friction  for  the  small  (10-1/2") 

model  flywheel.  Figure  18a  applies  to  forward-backward  motion  for  deflec¬ 
tion  about  the  upper  gimbal  axis  and  Figure  18b  to  forward-backward  motion 
about  the  lower  gimbal  axis.  Tlie  reason  for  these  separately  defined 
motion  planes  (90  degrees  apart)  is  that  the  gimbal  axes  do  not  intersect, 
thereby  causing  the  geometric  conditions  to  alternate  as  the  flywheel  whirls 
around.  This ,  along  with  the  dependence  on  more  system  parameters ,  makes  it 
seem  very  possible  to  overcome  the  instability  by  varying  the  system 
parameters . 
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Fig.  18a.  Internal  Friction  Force  For  Model  With  Bowed  Ropes 
For  Deflection  About  Upper  Gimbal  Axis 


The  stability  of  the  small  (10  1/2")  model  flywheel  has  been  observed 
to  be  slightly  enhanced  by  making  the  support  ropes  shorter.  To  illustrate 
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Fig.  18b.  Internal  Friction  Force  For  Model  With  Boived  Ropes 
For  Deflection  About  Lower  Gimbal  Axis 


how  the  inatheroatical  analysis  can  be  iised  to  investigate  this  sort  of  behav¬ 
ior,  the  coiiputer  program  was  used  to  predict  system  stability  for  support 
ropes  of  various  lengths .  Figure  19  shows  the  results  of  this  work  and  it 
can  be  seen  that  the  threshold  speed  of  stability  does  indeed  increase  ^len 
the  support  ropes  are  shortened.  There  is,  however,  a  practical  limit  to 
how  short  the  ropes  can  be  made  beca^Jse  synchronous  vibration  due  to 
unbalance  starts  to  increase  rapidly  as  the  ropes  are  made  shorter. 

III.  Computer  Solution 

The  differential  equations  of  motions  for  the  rigid  shaft  model  and  the 
bowed-out  ropes  model  are  given  in  Appendix  B.  These  equations  represent 
mathematical  relationships  between  all  the  forces,  external  and  internal, 
acting  on  the  system  and  the  resulting  motion  of  the  system  v/hich  at  this 
point  is  unlmoun.  For  each  model,  the  equations  form  a  linear,  eighth-order 
system  comprised  of  four,  linear,  second-order  differential  equations  whidi 
are  homogeneous  and  have  constant  coefficients.  These  equations  can  be 
solved  to  yield  the  general  solution  in  terms  of  some  unknown  constants  as 
follows. 

First,  the  four,  second-order  equations  can  be  converted  to  eight,  first- 
order  equations  by  the  definition  of  four  new  coordinates  which  are  equal  to 
the  velocities  of  the  four  existing  coordinates.  Substitution  of  these  four 
new  variables  into  the  second-order  equations  eliminates  the  second-order 
terms  and  reduces  the  equations  to  first-order.  This  procedure  yields  eight, 
first-order  equations  which  can  be  written  in  matrix  form  as 

f  =  [A]x  (1) 

where  x  is  the  state  vector  describing  the  configuration  of  the  system  and 
X  is  the  time  derivative  of  x.  The  matrix  A  is  the  square  matrix  made  up  of 
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Fig.  19.  Threshold  Speed  o£  Instability  for  10-1/2"  Flywheel 
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the  coefficients  of  tlie  four  original  equations.  The  solution  to  this  matrix 
equation  is  knoii/n  to  be  an  exponential  of  the  form 

X  =  C  e^^  ^2) 

where  C  is  a  vector  the  same  rajil\  as  x,  and  A  is  a  constant  knoun  as  an  eigen 
value.  In  general,  C  and  X  can  both  be  complex  having  real  and  imaginary  parts. 
Substituting  (3)  into  (ll  gives 

XCe^^  =  [A]Ce^^ 

or  [[A]  -  [I]A)  Ce^^  =  0  (5) 

For  this  equation  to  have  a  non- zero  solution,  it  is  necessary  that  the  matrix, 
defined  by  the  exyression  within  the  parentheses,  not  be  invertible.  ITiis  re¬ 
quires  that  its  determinant  be  equal  to  zero  ,  or 

|[A]  -  [I]x|=  0 

The  determinant  will  be  zero  only  for  special  values  of  X,  termed  eigeni allies. 

Tire  eigeni^alues  are  obtained  by  expanding  the  determinant  into  a  polinominal  in 
X,  and  then  solving  for  the  roots.  In  the  case  of  the  flmiheel,  the  resulting 
poliTiomial  will  be  eight-degree  and  will  thus  yield  eight  eigenvalues.  In  gen 
eral,  this  eight-degree  poliTiomial  is  impossible  to  solve  analytically  and  so 
the  computer  is  used  to  find  the  eigem^alues  using  a  program  outlined  in 

Appendix  C. 

Tne  eigenvalues  contain  some  of  the  most  important  information  about  the 
svstem.  For  each  particular  mode  of  motion  for  the  system,  there  is  an 
eigenvalue  that  pertains  to  it.  Hie  imaginan^  part  of  the  eigenvalue  is  the 
damped  natural  frequenc)-  for  the  particular  mode.  The  real  part  of  tlie  eigen¬ 
value  is  the  damping  exponent  and  it  determines  the  stability  of  motion  in  the 
mode.  If  this  damping  exponent  is  a  negative  number,  then  any  motion  in  this 
mode  will  decay  with  time  and  eventually  die  out  so  the  system  is  stable.  Should 
it  be  a  positive  number,  then  the  motion  will  grow  exponentially  with  time  and 
the  system  is  termed  unstable. 
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In  order  to  determine  what  shape  the  system  takes  on  when  there  is 
motion  in  a  particular  mode  (i.e.,  the  mode  shape),  it  is  necessary  to  first 
find  the  eigenvector  corresponding  to  the  mode  in  question.  The  eigenvector 
is  a  complex  vector  that  is  obtained  by  taking  an  eigenvalue  and  substituting 
it  into  Equation  (2)  and  solving  the  resulting  system  of  linear  algebraic 
equations  for  the  complex  constant  C  which  will  be  the  eigenvector.  The  actual 
mode  shapes  are  then  obtained  directly  from  the  eigenvectors  by  a  procedure 
discussed  in  Appendix  C. 

CONCEPTUAL  DESIGN  AND  COST  ANALYSIS  OF  A  50  KI\'h  FLEXIBLE  FLYIMEEL  ENERGY  STOR/\GE 
SYSTEM 

Tlie  perfomance  requirements,  which  form  the  basis  for  this  analysis, 
are  tliat  the  flvavheel  interfaces  with  a  small  solar  electric  energy  system 
for  residential  applications.  In  meeting  this  criteria,  the  fl>avheel  must  be 
sized  to  store  a  nominal  50  KWh  of  energy  and  the  motor/generator  must 
supply  10  KIV  peak  power. 

Definition  of  Components 

Rotor:  Tlie  rotor  is  comprised  of  a  flexible  flywheel  which  is  inherently 
self-balancing.  The  flyi\fheel  is  a  simple  hoop,  constructed  by  coiling  ropes, 
cables,  or  fibers  in  a  circle  of  the  desired  radius  until  the  design  mass  is 
achieved.  In  operation,  centrifugal  force  maintains  the  circular  configura¬ 
tion.  To  maintain  shape  integrity  when  not  spinning,  an  outer  sheath  surrounds 
the  cross  section  of  the  hoop.  The  hoop  is  suspended  by  flexible  ropes  which 
carr>^  only  the  weight  of  the  fl>avheel  and  provide  a  non-rigid  connection  be- 
tw'een  the  rotor  and  the  motor  generator. 

Bearings:  The  motor  bearings  provide  support  for  the  fl>’wheel  in  addition 
to  their  normal  functions  in  the  motor.  This  design  both  minimized  the  number 
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of  bearings  in  the  system  and  eliminated  the  necessity  for  a  clutch  between 
the  motor/generator  and  rotor.  These  bearings  must  be  capable  of  high 
speed,  low  drag  operation  in  order  to  minimize  system  losses. 

Vacuum  Vessel;  In  order  to  keep  air  drag  losses  down  to  acceptable 
levels,  the  rotor  assembly  must  operate  in  a  vacuum  vessel  capable  of  main¬ 
taining  a  vacuum  of  10  microns  or  less.  The  vacuum  vessel  performs  the 
the  additional  role  of  providing  support  assemblies  from  which  the  gimballed 

motor /generator  and  rotor  are  suspended. 

Vacuum  System:  In  order  to  maintain  the  vacuum  vessel  at  or  below 
10  microns,  two  vacuum  system  options  are  available  depending  on  the  leak 
rate  of  tlie vessel.  By  designing  the  vessel  with  a  low  leak  rate,  only 
occasional  evacuation  would  be  required.  This  option  emdsions  a  mobile 
vacuum  service  periodically  visiting  the  different  sites  and  performing  pump 
down  with  on-board  equipment  as  required.  The  second  option  is  to  have  a 
vacuum  system  pemanently  attached  to  the  vacuum  vessel  providing  continuous 
evacuation.  This  option  allows  higher  leak  rates  at  the  cost  of  energy  to 
run  the  pumps. 

Electronics:  The  electronics  provide  the  interface  between  the  motor/ 
generator  and  input-output  power  requirements .  On  input  the  electronics  must 
perform  the  task  of  maximum  power  tracking  to  allow  efficient  transfer  of 
power  from  the  solar  array  to  the  motor.  On  output  the  electronics  system 
must  accept  the  electrical  energy  from  the  generator  and  cariy  out  the  re¬ 
quired  power  conditioning  (DC  to  AC,  voltage  changes,  etc.)  to  allow  efficient 
use  at  the  load. 

Safety  Devices:  The  flexible  flywheel's  gradual  failure  mode,  unraveling 
allows  for  a  minimum  of  safety  devices.  Sensors  in  the  vacuum  vessel  can  moni¬ 
tor  for  early  warning  signs  such  as  hoop  growth  or  whiplash  and  shut  down  the 
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flywheel  system  before  final  failure.  Even  if  failure  should  occur,  the  rotor 
does  not  rapidly  disintegrate  due  to  its  self-balancing  characteristics.  As 
such,  safety  devices  need  not  be  as  rigorous  and  extensive  as  with  solid 
rotors . 

Oil  cooling  system:  To  provide  cooling  for  the  motor/ generator ,  trans¬ 
former  oil  is  circulated  around  the  motor  windings,  llie  system  consists  of 
a  heat  exchanger,  oil  pump  and  cooling  lines  to  the  motor  generator.  The 
cooling  lines  inside  the  vacuum  vessel  must  be  flexible  to  allow  the  gimbal 
action  of  the  motor/generator. 

Elements  Defined  By  Technological  Requirements 

Three  components  of  the  flyivheel  system  are  defined  by  technological 
constraints.  These  elements  are  the  rotor,  the  moto r/ generator ,  and  the 
bearings.  These  requirements  and  their  relationship  to  each  component  are 
given  in  Table  3. 

A  unique  feature  of  the  "flexible  flyv\'heel"  is  the  construction  of  the 
rotor  from  flexible  fibers  (such  as  synthetic  rope)  with  no  bonding  agents. 

The  advantages  of  this  configuration  are: 

1.  Self  balancing-highly  flexible  rotor  operates  at  supercritical  speeds. 

2.  Simple  construction. 

3.  High  strength  fibers  in  pure  tension  with  no  bonding  materials  to 
create  mismatch  in  elasticity  or  strength. 

4.  Gradual  failure  mode,  with  early  warning  by  whiplash  or  hoop  growth. 

Figure  20  shows  tlie  basic  design  of  the  rotor  hoop  as  well  as  some  practical 

limits  on  the  radii  involved.  The  height  (L)  of  the  hoop  is  varied  with  each 
material  in  order  to  obtain  the  rotor's  required  mass.  The  practical  limit  on 
L  is  probably  about  four  feet.  A  large  L,  however,  introduces  an  additional 
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TABLE  3.  Satisfaction  of  System  Requirements 


Technological  Requirement 

Component  affected/ 

Solution  for  Problem 

50  KlVh  storage  capacity 

Self-balancing 

Gradual  failure  mode 

Failure  early  warning 

Hoop  height  less  than  4'  for 
stability 

Flexible  hoop  rotor 

High  drag  in  rotary  vacuum  feed¬ 
through 

Gimballing  action  required  at 
motor  for  stability 

Motor/generator  in  vaaium 

Need  low  support  stiffness 

Need  asymmetric  support  stiffness 

Need  support  damping  and  mass 

Gimbal  motor/generator  on 
nonintersecting  axis 

Low  drag 

High  speed  operation 

Bearings 
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Fig.  20.  Rotor  Hoop  Dimensions 
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dynamic  whirl  mode  with  its  potential  instability.  Therefore,  the  height  of 
the  rotor  should  be  kept  within  reasonable  limits. 

Possible  materials  for  the  construction  of  the  rotor  include  Nylon,  Dacron, 
high  strength  carbon  steel  cable,  and  Kevlar- 29  fiber.  A  comparison  of  rotor 
hoops  constructed  of  these  materials  can  be  found  in  Table  4.  Based  on  keeping 
the  rotor  height  to  a  minimum,  one  can  quickly  eliminate  Nylon  and  Dacron  from 
further  consideration.  Carbon  steel  cable  and  Kevlar  fiber  have  acceptable 
hoop  height  requirements  although  Kevlar  represents  almost  70%  reduction  of  hoop 
height  when  compared  with  steel.  However,  the  loxv  energy  storage  density  of 
carbon  steel  cable  makes  it  unattractive  as  a  candidate  for  construction  of  a 
flexible  rotor. 

A  hoop  constructed  of  Kevlar- 29  capable  of  storing  50  KlVh  would  weigh 
about  810  pounds  and  operate  with  a  maximum  rpm  of  14,600.  The  resulting 
rotor  would  have  a  3'  major  radius  and  would  be  a  nominal  8”  thick  by  about 
10.4  inches  in  height.  Based  on  material  costs,  the  1500-Den  Kevlar  29  fiber 
would  be  tlie  best  choice  among  the  various  sizes  available. 

It  has  been  shown  that  in  order  to  overcome  the  subsynchronous  whirl  in¬ 
stability  of  the  flexible  flywheel,  the  motor/ generator  must  be  gimballed  on 
nonintersecting  axes  to  provide  tlie  inportant  parameters  for  whirl  stability, 
namely,  low  support  stiffness,  stiffness  asymmetry,  and  bearing  support  mass 
(the  motor  itself) .  In  order  to  place  the  motor/ generator  on  the  outside  of 
the  vacuum  vessel,  a  vacuum  feedthrough  with  three  degrees  of  rotational  free¬ 
dom,  mvo  for  the  gimballing  action  plus  one  for  rotation  of  the  rotor,  would  be 
required.  Such  a  feedthrough  would  be  costly  to  design  and  build  and  could 
adversely  affect  the  whirl  stability  parameters.  Tlierefore,  as  a  result  of 
the  gimballing  requirements,  the  motor/ generator  must  be  placed  inside  the  vacuum 
vessel.  Since  a  number  of  motor/ generators  have  been  designed  within  the  last 
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TABLE  4.  Calculated  Values  for  50  KWh  6- foot  Diameter  Flywheel  (Fig.  1) 
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Safety  factor  = 


decade  for  high  speed  operation  in  a  vacuum,  this  requirement  should  present 
no  major  problem. 

Additional  requirements  placed  on  the  motor/generator  are  10  KtV  peak 
power  output  and  constant  power  operation  over  at  least  a  2:1  speed  range. 

If  it  were  not  for  the  requirement  of  operation  in  a  vacuum,  hydrody¬ 
namic  (fluid  film)  bearings  would  be  ideal  for  energy  storage  flywheel  appli¬ 
cations,  However,  the  lubricating  fluid  required  would  vaporize  to  some 
extent  and  necessitate  greater  pumpdown  requirement.  Any  seal  which  might 
be  iLsed  to  overcome  this  problem  would  add  unacceptable  amounts  of  power  loss 

thereby  reducing  the  useful  energy  storage  period. 

A  major  disadvantage  of  ball  bearings  for  flywheel  applications  is  the 
cage  and  retainer  instabilities  which  appear  at  high  speeds.  Research  re¬ 
ported  by  Eisenhaure  and  Kingsbury  (Reference  4)  has  shown  that  this  dis¬ 
advantage  can  be  overcome  by  using  retainerless  bearings.  Thus,  the  retainer¬ 
less  precision  ball  bearing  appears  to  be  well  suited  to  the  energy  storage 
fl>nvheel  system.  However,  these  bearings  also  require  a  light  fihii  of  lubri¬ 
cant  for  operation. 

Magnetic  bearings  offer  the  advantage  of  needing  no  lubricant  and  have  been 
shown  to  have  ver)^  low  drag  and  long  life.  Power  for  the  bearing  can  be  supplied 
by  a’uxiliary  windings  on  the  motor  generator  allowing  failsafe  spin  down  opera¬ 
tion.  Mechanical  touchdown  bearings  would  only  be  needed  for  cold  start/stop 
situations.  Although  magnetic  bearings  are  more  costly  than  other  alternatives, 
the  cost  would  more  than  be  offset  over  the  life  of  the  system  in  decreased 
maintenance  costs. 

Elements  Subject  to  Economic  Trade-offs 

Miile  the  need  for  operating  a  higher  performance  energy  storage  flywheel 
in  a  vacuum  is  not  subject  to  economic  trade-offs,  the  method  of  providing  and 
maintaining  the  vacuum  is.  These  trade-offs  affect  the  vacuum  vessel,  vacuum 
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system  aiid  as  mentioned  previously,  the  use  of  lubricated  or  oil  free  bearings. 

Hie  vacuum  vessel  perfonns  a  dual  role:  1)  to  provide  a  low  drag  environ¬ 
ment  for  the  rotor;  and  2)  to  provide  a  support  structure  from  which  to  suspend 
the  motor/generator  and  rotor  assembly.  In  addition,  if  the  system  is  installed 
above  groimd,  the  vacuum  vessel  must  provide  safety  containment  should  a  cata¬ 
strophic  failure  occur.  Two  major  options  exist  for  the  constiaiction  vessel, 
steel  or  reinforced  concrete. 

Hie  steel  vessel  is  designed  to  provide  ease  in  final  assembly  and  install¬ 
ation  of  the  system.  Hie  vessel  is  comprised  of  two  sections.  The  lower  sec¬ 

tion  consists  of  a  steel  cylinder  having  semihemi spheric  bottom  and  a  ring 
flange  on  the  top  to  mate  with  the  upper  section  ivhich  is  also  semihemispherical . 
It  contains  all  the  necessar)^  feedthroughs  and  support  structures  for  the  motor/ 
generator  and  rotor  assembly.  Also  provided  on  the  upper  section  is  the  pumpdown 
port.  Hiis  arrangement  allows  the  entire  system  to  be  assembled  on  the  upper 
section,  then  lowered  into  the  lower  section  of  the  vacuum  vessel.  Hie  two  sec¬ 
tions  are  then  welded  together  at  the  mating  flanges.  Hiis  design  also  has  the 

advantage,  should  the  motor/generator  or  rotor  fail,  only  the  upper  section  need 
be  removed,  by  cutting  the  weld,  for  repair  or  replacement. 

Hie  steel  vessel  offers  the  advantage  of  low  outgassing  and  very  low  poro¬ 
sity,  thus  decreasing  the  vacuum  maintenance  requirements.  Ulien  installed  below 
ground  the  steel  cylinder  need  only  be  1/4”  thick  with  the  surrounding  earth 
providing  failure  containment.  However,  when  installed  above  ground,  the  thick¬ 
ness  must  be  increased  to  5/8”  to  provide  sufficient  containment  during  failure. 

The  concrete  vessel  uses  the  same  upper  head  assembly  with  its  associated 
advantages.  Hie  lower  section,  however,  is  constructed  of  2”  thick  reinforced 
concrete  with  a  thin  outer  shell  of  steel.  Hie  outer  shell  provides  the  mating 
flange  for  welding  to  the  upper  head.  Hie  shell  is  also  needed  because  of  the 
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high  porosity  in  the  concrete,  which  would  allow  air  to  diffuse  into  the  vacuum 
chamber.  This  shell  need  only  be  about  1/32"  thick  with  the  concrete  providing 
all  the  needed  structural  strength. 

The  main  advantage  of  tiie  concrete  vessel  is  realized  in  above  ground 
installation  where  the  concrete  provides  the  needed  failure  containment.  Un¬ 
fortunately,  concrete  outgasses  in  a  vacuum. raising  pumping  requirements  to 
the  point  v\fhere  continuous  on-site  pumping  is  required. 

In  order  to  keep  windage  drag  down  to  acceptable  limits,  a  vacutun  of  better 
than  10  microns  (lO'^  torr)  is  required.  To  maintain  this  vacuum,  the  pumping 
system  can  be  either  on-site  (continuous)  or  brought  to  the  site  as  needed 
periodically.  Continuous  on-site  pumping  extracts  a  price  of  higher  initial 
capital  outlay  for  the  system.  In  addition,  energy  must  be  supplied  to  the 
punp.ing  system  for  operation  thereby  lowering  the  overall  system  efficiency. 
Tlierefore,  the  on-site  system  should  only  be  used  where  required. as  in  the  con¬ 
crete  vessel  option.  The  on-site  vacuum  system  would  consist  of  a  mechanical 
roughing  pump  just  large  enough  to  overcome  the  leaks  and  outgassing  at  the  10 
micron  level.  Initial  pump  down  would  be  done  at  the  time  of  installation  by 
a  unit  brought  to  the  site  for  that  purpose. 

In  the  periodic  pumping  sdieme,  the  vacuum  system  would  be  brought  to  the 
site  twice  a  year.  A  nationwide  netitfork  of  companies  could  provide  this  service 
as  is  presently  done  with  many  other  services  nowadays.  The  service  equipment, 
consisting  of  a  truck  containing  a  diffusion-pump  based  vacuum  system,  would  be 
brought  to  the  site  and  hooked  up  to  tlie  vacuum  port.  The  system  would  then  be 
pumped  dovn  to  lO’^  torr  and  valved  off.  Hhen  the  pressure  reached  10  ^  torr, 
the  process  needs  to  be  repeated. 

Tlie  advantages  of  a  permanent  magnet  (PM)  DC  mo  tor/ generator  design  for 
high-speed  applications  are  high  efficiency  and  relatively  simple  power  condition¬ 
ing  electronics.  The  next  best  candidates  are  the  functionally  similar  Lundell 
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and  Nadyne  homopolar  inductor  designs  with  the  field  supplied  by  a  stationary 
coil,  and  without  brushes  or  slip  rings.  However,  in  addition  to  field  coil 
paver  dissipation,  the  latter  designs  are  heavier  and  bigger  than  their  PM 
counterparts  and  may  therefore  produce  side  loading  forces  which  are  a  problem 
for  magnetic  bearing  systems.  The  advantage  of  a  controllable  field  over  a 
PM  field  is  often  offset  by  the  need  for  more  complex  and  expensive  power  con¬ 
ditioning  electronics.  The  primary  disadvantage  of  PM  units  for  energy  storage 
devices  is  the  electromagnetic  drag  due  to  stator  eddy  currents. 

A  brushless -DC  PM  motor/ generator,  under  development  by  A.  R.  Millner  of 
MIT/Lincoln  Laboratory  (Reference  3) ,  avoids  the  eddy-current  problem  by  using 
a  truly  ironless  armature.  The  only  conducting  material  in  the  stator  is  the 
copper,  and  this  is  made  of  small  conductors  individually  insulated  for  minimum 
eddy  current  losses.  Figure  21  shows  the  layout  of  the  Lincoln  Laboratory  motor/ 
generator  with  magnetic  bearing.  Tliis  motor/generator  design  appears  to  be  a 
good  match  for  the  dual  tasks  of  solar  energy  storage  and  conversion.  Used  in 
conjunction  with  the  PM  motor/generator,  a  naturally  coiinnutated  DC  input  bridge 
can  closely  track  the  maximum  power  point  of  the  energy  source,  thereby  eliminat¬ 
ing  mismatches  that  lower  system  efficiency.  Furthermore,  using  a  cy do converter, 
a  good  quality  60  Hz  AC  output  for  the  load  can  be  produced. 

Even  through  MIT's  motor/generator  design  has  high  efficiency  (on  the  order 
of  90%)  it  is  felt  that  more  cooling  will  be  needed  than  can  be  provided  by  con¬ 
duction  through  the  frame,  gimbal  and  support  structure.  This  is  provided  by 
"canning"  the  stator  armature  and  flowing  transformer  oil  around  the  windings. 

Tlieoil  cooling  system  consists  of  a  heat  exchanger,  pump  and  flexible  oil 
lines  to  the  motor  (so  as  not  to  inhibit  the  gimballing  action) .  By  assuming 
1  KIV  of  heat  must  be  removed  (based  on  90%  efficiency  10  KIV  input  output)  ,  a 
flow  rate  of  .25  gal/min  would  be  sufficient  to  limit  the  temperature  rise  in 
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the  oil  to  less  than  20°C.  A  small  oil  pimp  consuming  less  than  200  watts 
would  be  sufficient  to  provide  the  needed  flow  rate. 

llie  heat  exchanger  offers  two  options,  air  cooled  or  water  cooled.  The 
air  cooled  option  offers  the  advantage  of  being  a  passive  heat  exchanger,  the 
cooling  being  provided  by  natural  convection.  Unfortunately,  this  necessitates 
a  much  larger  surface  area  to  dissipate  tlic  heat.  Hie  water  cooled  unit,  while 
smaller  in  size,  requires  a  ready  source  of  water  such  as  a  lake  or  stream.  If 
the  water  is  available,  this  option,  by  far,  would  be  the  dioice;  however,  if 
the  water  must  be  supplied  by  a  well,  the  energy  expended  by  pumping  the  water 
would  be  unacceptable. 

The  question  arises,  sliould  the  energy  storage  system  provide  the  output 
in  the  form  of  AC  or  DC  power?  The  AC  option  offers  the  advantage  of  not  having 
to  make  any  change  at  the  user  end  of  the  system.  By  providing  60  Ilz  220V  power 
to  residential  users,  the  system  could  be  backfitted  easily  to  present  home 
systems.  Also,  when  backup  power  from  a  utility  was  needed,  a  simple  switch¬ 
over  could  be  accomplished.  Tlie  disadvantage  would  be  the  cost  of  the  cyclo¬ 
converter  needed  to  convert  the  output  of  the  flywheel  to  AC. 

On  the  other  hand,  DC  output  could  provide  most  of  the  energy  needs  in 
the  average  home.  All  heat  producing  products,  such  as  stoves,  water  heaters, 
furnaces,  etc.,  could  run  on  DC  power.  Inverters  for  units  requiring  AC  power 
such  as  motors  and  compressors  or  electronic  devices  like  T.V.  could  be  provided 
individually.  However,  this  creates  the  problems  of  a  dual  electrical  system, 
one  for  DC  operating  units  and  the  other  for  the  AC  units.  Also,  should  backup 
power  be  needed  from  a  utility,  it  would  have  to  be  rectified  and  conditioned 
to  run  in  the  DC  electrical  units. 

Let’s  now  consider  two  possible  systems  based  on  the  previous  design  require¬ 
ments.  System  A  would  be  an  above  ground  installation  and  System  B  would  be 
below  ground. 
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System  A  would  consist  o£  a  concrete  vessel  sitting  above  ground.  This 
would  be  a  major  advantage  in  rocky  areas  where  below  ground  installation 
would  be  prohibitively  expensive.  As  mentioned  previously,  this  system 
would  require  an  on-site  continuous  vacuum  system.  For  comparative  purposes, 
this  system  will  utilize  the  DC  electrical  output  option. 

System  B  is  a  below  ground  steel  vessel  unit.  The  system  utilizes  a 
periodic  pumpdown  vacuum  scheme  and  AC  electrical  output.  Table  5  gives  the 
comparable  system  costs  of  the  options  for  production  model  quantities  of 
10,000  units  per  year. 

The  motor/generator,  magnetic  bearings,  and  associated  electronics 
costs  were  obtained  from  Reference  3.  These  costs,  based  on  1979  technology , 
are:  1)  the  motor/generator,  $1000  @  $100/KW;  2)  the  magnetic  bearings, 

$3750  @  $75/KWh;  3)  the  electronics,  $1600.  The  costs  of  the  two  systems 
(1979  $)  are  coirparable;  System  A  $24,241  ($485/KWh)  and  System  B,  $24,160 
($483/KWh) . 

To  get  a  better  comparison  of  the  two  systems,  a  cost  analysis  based  on 
the  Life -Cycle  Cost  Methodology  of  Reference  5  was  performed.  Input  for  the 
analysis  is  shown  in  Table  6. 

The  methodology  computes  the  residential  total  present  value  (1985  year 
of  operation  selected)  and  total  annual  costs  including  all  the  parameters 
listed  in  Table  6.  The  cost  analysis  reveals  the  following  values: 

System  A  System  E 

Total  Present  Value  (1979  $)  $44,642  $36,437 

Total  Annual  Cost  $  5,244  $  4,280 

The  cost  difference  between  the  two  systems  primarily  arises  from  the  higher 
initial  and  annual  costs  of  a  continuous  on-site  vacuum  system. 
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T.ABLE  5.  Cost  Breakdown  (1979  Dollars)  o£  System  A  and  B  (50  Kl\li)  for 
Production-^'bdel  Quantities  of  10,000  units/year 


Rotor 

Metor/ Generator 

Bearings 

Gimbal 

Support  Structure 
Input  Electronics 
Output  Electronics 
Inverters 

Installation  and  Assembly 
Oil  Cooler 
Vacuum  Vessel 
Vacuum  System 
Backup  Power  System 

Total 

One  Year  Maintenance 


System  A* 

System  B* 

$  6,885 

$  6,885 

1,000 

1,000 

3,750 

3,750 

2,125 

2,125 

1,400 

1,400 

800 

800 

50 

800 

774 

0 

1,500 

2,000 

200 

200 

3,580 

4,850 

1,500 

300** 

650 

50 

$24,241 

$24,160 

1,000*** 

500** 

*  Based  on  a  2,000  sq.  ft.  home  with  air  conditioner,  refrigerator,  freezer, 
and  electric  furnace. 

**  Vacuum  system  costs  including  vacuum  valve  and  system  maintenance. 

***  Includes  cost  of  operating  vacuum  pump. 
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TABLE  6.  Input  for  Cost  Analysis  by  Life-Cycle  Cost  Methodolog)^** 


System  A 

System  B 

Discount  rate  (d)* 

.10 

.10 

Percent  down  payment  (D)* 

.20 

DO 

O 

Energ)^  escalation  rate  (e)* 

.o:' 

.07 

Backup  energ>"  cost  (F) 

- 

- 

General  inflation  rate  (g)* 

.05 

.05 

Interest  rate  (i)* 

.115 

.115 

Investment  tax  credit  (I)* 

0 

0 

Initial  system  cost  (ICp) 

$24,241 

$24,160 

Period  of  analysis  (N)* 

20  yrs 

20  yrs 

Accounting  lifetime  (N^)* 

- 

Borrowing  period 

30  yrs 

30  yrs 

Miscellaneous  (maintenance)  cost  (OM) 

$  1,000 

$  500 

Property  taxes  (P)* 

0 

0 

Salvage  value  (S\0  * 

0 

0 

Income  tax  rate  (t)* 

.30 

.30 

Year  of  operation  (yo) 

1985 

1985 

Price  year  (>Tj) 

1979 

1979 

*  Values  from  ”M  Analysis  of  the  Current  Economic  Feasibility  of  Solar 

Water  and  Space  Heating,"  Office  of  the  Assistant  Secretan' 
for  Conservation  and  Solar  Applications,  DOE/CS-0023.  Jan- 
uar>'  1978 

**  A.  M.  Perino,  "A  Methodology  for  Determining  the  Economic  Feasibility  of 
Residential  or  Commercial  Solar  Energy  Systems,"  Sandia 
Laboratories,  SAND  78-0931,  January  1979. 
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Suinmary  of  Best  Option 

Hie  conceptual  design  of  a  50  KWh  Flexible  Flywheel  Energy  Storage  System 
is  shown  in  Figure  22.  The  system  is  designed  for  below  ground  installation 
although  above  ground  or  partially  below  ground  installation  would  represent 
only  minor  changes  and  not  greatly  affect  the  cost  of  the  system.  The  system 
stores  a  nominal  50  KlVh  with  a  peak  output  of  10  Kb'  AC. 

The  major  components  of  the  system  are:  1}  the  rotor;  2)  the  gimbal 
mounted  motor/generator,  magnetic  bearing  and  electronics;  5)  the  oil  cooling 
system;  and  4)  the  vacuum  vessel. 

Tlie  rotor  is  a  simple  hoop  constructed  by  coiling  Kevlar  fibers  into  a 
circle  6'  in  diameter  with  a  minor  radius  of  4".  To  maintain  shape  integrity 
when  not  spinning,  an  outer  rope  sheath  surrounds  the  cross-section  of  the 
hoop.  The  overall  hoop  assembly  weighs  850  pounds  and  stores  a  nominal  50  Kl\h 
at  14,000  rpm.  Maximum  rpm  is  14,600  with  a  safety  factor  of  2.  Four  flexible 
ropes  attach  the  hoop  to  the  motor/ generator  shaft. 

Hie  oil  cooling  system  consists  of  a  radiator,  oil  pump  and  flexible  lines 
to  the  mo tor/ generator.  The  radiator  is  an  oil-to-air  cooling  unit  similar  to 
those  found  in  large  cars  and  trucks.  A  small  oil  pump  is  included  to  provide 
flow  in  the  systems.  Because  the  pressures  involved  would  be  low,  the  pump 
ivould  only  have  to  overcome  small  head  losses  and  only  require  50  to  100  watts 
of  power.  The  flexible  lines  are  needed  to  allow  free  movement  of  the  motor 
generator  in  its  gimbal  mounts.  The  system  uses  insulating  transformer  oil 
to  eliminate  the  possibility  of  shorting  out  the  armature  stator. 

The  vacuum  vessel  is  designed  to  provide  ease  in  final  assembly  and  in¬ 
stallation  of  the  system.  The  vessel  is  comprised  of  two  sections.  The  lower 
section  consists  of  a  1/4”  thick  steel  cylinder  with  a  semihemispherical  bottom. 
The  upper  section,  also  semihemispheric,  contains  all  necessary  feedthroughs  for 
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the  electric  oil  lines  and  the  support  structure  for  the  motor  generator 
and  rotor  assembly.  Also  provided  is  a  pumpdovm  port  and  vacuum  valve. 

This  arrangement  allows  the  entire  system  to  be  assembled  on  the  upper  head 
tlien  lowered  into  the  rest  of  the  vacuum  vessel.  The  two  sections  are  then 
welded  together  at  the  mating  flanges  either  at  the  manufacturing  plant  or 
at  the  installation  site.  This  design  also  has  the  advantage,  should  the 
motor -generator  or  rotor  fail,  that  only  tlie  upper  section  need  be  removed, 
by  cutting  the  weld,  for  unit  repair  or  replacement  instead  of  removing  the 
entire  vacuum  vessel  from  its  below  ground  installation. 

A  motor/generator  similar  to  the  one  being  developed  at  MIT/Lincoln 
Laboraton^  (A.  R.  Millner)  will  be  used.  A  significant  departure  from 
Millner's  design  is  to  "can"  the  stator  armature  to  enable  cooling  of  the 
windings  with  transformer  oil.  The  MIT  motor/generator  is  being  developed 
under  the  sponsorship  of  the  U.  S.  Department  of  Energy  (DOE) ,  specifically 
for  flwheel  applications  interfaced  with  solar  photovoltaic  systems. 

Tie  selected  motor/generator  is  of  a  brushless -DC,  permanent  magnet 
design  offering  the  advantages  of  high-efficiency  and  relatively  simple  power 
conditioning  electronics  for  high  speed  applications.  Oil  cooling  of  the 
stator  armature  was  added  to  the  design  because  it  is  believed  that  conduction 
through  the  frame  and  gimbal  support  structure  would  provide  insufficient  cool¬ 
ing.  Futher  detail  on  the  motor/generator  as  well  as  the  bearings  and  elec¬ 
tronics  can  be  found  in  Reference  3. 

The  magnetic  bearings  were  chosen  over  ball  bearings  to  avoid  the  neces¬ 
sity  of  bearing  lubrication.  In  the  case  of  ball  bearings,  the  oil  film  would 
constantly  be  boiling  off  into  the  vacuum  chambers,  necessitating  constant  oil 
replacement  and  more  frequent  vacuum  vessel  pumpdoi'm,  both  resulting  in  higher 
maintenance  costs.  Tliese  maintenance  cost  increases  would  more  than  offset. 
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over  the  life  of  the  system,  the  higher  cost  of  the  magnetic  bearings  (based 
on  the  life-cycle  costing  methodology  used  in  the  cost  analysis) .  The  bear¬ 
ings  can  be  pov\fered  by  auxilliary  windings  on  the  motor,  allowing  fail¬ 
safe  spin-douTi  operation.  ^fecIlanical  touchdown  bearings  are  included  for 
cold  start/stop  situations. 

The  electronics  interface  between  the  DC  solar  array  and  the  PM  motor/ 
generator  provides  liigh  quality  60  Hz  AC  for  residential  power  consumption. 

All  solar/electric  power  would  go  into  the  motor/generator  to  spin  up  the 
rotor.  Tlie  motor/generator  would  be  electronically  controlled  to  serve  as 
a  maximum  power  tracker,  Hus  is  important  because  the  varying  electrical 
output  of  solar  arrays  is  generally  mismatched  to  the  characteristics  of  the 
storage  system  and  load,  causing  inefficient  operation.  The  generator  out¬ 
put  is  fed  into  a  silicon  controlled  rectifier  (SCR)  cycloconverter  which, 
with  a  transformer,  provides  220V/60  Hz  output  to  the  load. 

Table  7  gives  a  comparative  cost  breakdown,  in  1979  dollars,  of  a  10  Kl\'h 
prototype  fludieel  system  and  the  50  KHTi  system  B  production  model  of  Table  6. 
Major  cost  items  are:  1)  the  rotor;  2)  the  motor/generator;  3)  the  magnetic 
bearings;  4)  the  electronics  and  5)  the  vacuum  vessel.  The  rotor  is  constructed 
of  1500-Den  Kevlar  at  a  cost  of  $8.00  per  pound.  This  results  in  a  cost  of 
$6800  (850  lbs.)  for  the  50  KhTi  (production  model)  rotor  and  $1296  (162  lbs.) 
for  the  10  KTdi  prototype.  Rotor  constuiction  costs  are  estimated  to  be  $0.10 
per  pound  or  $85  for  the  production  model  and  $1,00  per  pound  or  $162  for  the 
prototype.  The  total  rotor  costs  are  then  $6885  for  the  50  Kldi  rotor  and  $1458 
for  the  10  KlVh  prototype  rotor.  It  should  be  noted  that  a  large  savings  can 
be  realized  if  the  cost  of  Kevlar  decreases  in  the  future  ($850  for  each  dollar 
reduction  in  the  price  per  pound.) 
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TABLE  7.  Cost  Breakdown 

(1979  Dollars)  for  Prototype 

and  Production  Model 

10  (Prototype) 

50  KlVh 

(10,000  unit/year) 

Rotor 

$  1,458 

$  6,885 

Motor/ Generator 

1,000** 

1,000* 

Electronics 

1,600** 

1,600* 

Bearings 

1,250** 

3,750* 

Vaaium  Vessel 

5,150 

5,150 

Gimbal 

850 

2,125 

Support  Structure 

500 

1,400 

Installation  5  Assembly 

1,000 

2,000 

Oil  Cooler 

100 

200 

Backup  Power  Switchover 

- 

50 

$12,680 

$24,160 

Prototype  Engineering 
Services  $1,500.00 

Maintenance  $500/yr 
based  on  2  pumpdowns 

Testing  $1,500.00 

per  year 

*  A.  R,  Millner,  "A  Flywheel  Energy  Storage  and  Conversion  System  for  Solar 
Photovoltaic  Applications,”  ASME,  79-501-1,  March  1979. 


**  Prices  higher  per  KlVh  due  to  one  of  a  kind  prototype. 
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The  costs  per  KWh  for  the  50  KWh  production  model  are  significantly 
lower  than  for  the  10  KWh  prototype  due  to  the  dual  benefits  of  scale  and 
quantity  production-- $20  vs.  $100  for  the  motor/generator,  $75  vs.  $125  for 
the  bearings  and  $32  vs.  $160  for  the  electronics.  Again,  these  areas  are 
subject  to  cost  reductions  with  advances  in  technological  development. 

The  cost  of  the  vacuum  vessel  are  determined  based  on  1979  costs  of 
steel  vessels  of  similar  dimensions.  This  cost  is  $5,150  and  quite  subject 
to  fluctuations  in  the  price  of  steel. 
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CONCLUSIONS 


The  principal  conclusions  from  this  research  can  be  summarized  as  follows : 

(1)  Subsynchronous  whirl  instability  at  supercritical  speeds  due  to 
internal  friction  is  strongly  suppressed  by  a  gimballed  motor/ 
generator  support  system.  (This  finding  is  applicable  to  more 
conventional  flywheels  as  well.) 

(2)  A  flexible  flywheel  made  from  synthetic  rope  or  fiber  with  no 
bonding  agent  is  easy  to  construct  and  can  be  made  self¬ 
balancing  by  a  properly  designed  support  system. 

(3)  The  flexible  flywheel  system  offers  a  relatively  inexpensive  and 
efficient  tool  for  investigating  flywheel  rotor  dynamics  in 
general ,  since : 

(a)  failures  are  not  catastrophic 

(b)  critical  speeds  are  so  low  that  whirl  modes  can  be 
visually  observed,  and 

(c)  support  rope  configurations  can  be  varied  to  approxi¬ 
mate  the  shaft  constraints  of  conventional  flywheels. 

(4)  The  simplicity  of  the  flexible  flywheel  construction  promotes 
economical  production  using  conventional  manufacturing  processes 
and  equipment. 
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APPENDIX  A 


Flexible  Fl)^eel  Sizing  and  Performance  Calculation 
1.  On  the  Breaking  Speed  and  1-laximum  Storable  Energy  in  a  Stranded  Hoop  Flywheel 


The  maximum  energy  E  which  can  be  stored  in  a  stranded  hoop  flywheel  without 
breaking  it  is  a  function  only  of  the  mean  hoop  radius  R  and  the  total  tensile 
proof  load  P  of  the  hoop  fibers  (cumulative  strength) . 

The  maximum  speed  corresponding  to  E,  however,  is  also  a  function  of 
the  lineal  weight  per  unit  length  w  of  the  hoop  material. 

The  design  equation  for  maximum  storage  energy  is: 

E  =  ttRP  ft- lb. 
or 


E  = 


ttRP _ 

(2.654) (10^) 


KWh 


where  R  =  ft. 
P  =  lb. 

The  maximum  speed  is  given  by 


where  g  =  32.2  ft/sec^ 
w  =  Ib/ft. 

It  can  be  seen  that,  although  the  material  weight  w  does  not  affect  the 
maximum  energy  E,  a  practical  advantage  of  a  heavy  material  is  that  it  develops 
its  maximum  energy  at  a  lower  speed,  thus  reducing  windage  and  friction  losses. 
Also,  currently  available  electrical  machinery  is  speed  limited. 

The  dimensions  of  a  flywheel  for  home  energy  storage  are  limited  by  space 
available  for  the  unit.  A  practical  maximum  radius  (R)  seems  to  be  about  3  feet. 
If  the  hoop  cross  section  radius  is  more  than  about  4  inches,  the  flywheel  is  no 
longer  a  ''hoop'',  the  equations  above  do  not  apply,  and  the  safety  advantages  of 
individual  independent  strands  are  probably  lost  by  some  required  bonding  or  ty¬ 
ing  of  the  fibers.  Thus  the  cross  sectional  area  of  a  hoop  with  circular  cross 
section  is  limited  to  about  16  frin^. 
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For  a  flwheel  of  3  feet  radius  and  given  material,  the  only  available 
wav  to  increase  the  cumulative  proof  load  P  (and  E)  is  to  increase  the  ax^l 

len<7t]i  L  o£  the  hoop  (thus  making  the  flywheel  a  long  hollow  cylinder).  Th 
practical  limit  on  L  is  probably  about  4  feet.  A  long  L,  however,  intro^ces 
another  dynamic  whirl  mode  into  the  system, _with  its  potential  instabi  i  y. 
Therefore  L  should  be  kept  as  short  as  possible. 

2.  Calculation  of  Size  and  Performance  in  Terms  of  Hoop  Radius  R 

In  sizing  a  flexible  fl>^heel,  it  is  helpful  to  have  performance  and  design 
parameters  expressed  in  terms  of  tlie  hoop  radius  R. 

The  following  calculations  give  the  maximum  energy  stored  E ,  corresponyng 
speed  Ng,  hoop  weight  W,  energy  density,  cost,  cost/Kl%,  and  bearing  power  loss 

in  terms  of  the  hoop  radius  R.  Also  given  are  the  y solute  values  for  a  2  ft. 
radius  hoop,  as  an  example  for  comparison  of  the  effect  of  materia.  . 

Tlie  calailations  are  made  for  DACRON^  1/4"  rope,  KEVLAR  1/4"  rope  (hand 

construction) ,  and  KEITAR^  -  49-1429  Denier  fiber  (machine  construction) .  All 
calculations  contain  a  factor  of  safety  -  2.0. 


A,  DACRON^  hheel  (1/4"  (t)rope)  ,  F.S.  =  2 

For  8"  (j)  cross  section  (W  =  .022  Ib/ft): 


50.26  in' 


P  ^  ,  =  1800  lb. 

strand 

M  -  50.26  . 

-^strands  .049  ° 


density g^o  strands,  xvhere  .049  is  c/s  area  per  strand 
factor  /  - 


P  =  (820) (1800)  =  1.476  x  10  lb, 


KTT  -y 

_ i ^  =  (.8735)R  KlVh,  where  P/2  reflects  a  factor  of  safety  of  2 

(2.654) (10°) 


M  -  60 

^E  27rR 


IV  =  lirRN  1  w  =  r  2tt  )  (820)  (.021  Ib/ft)  R  -  113  R  lb 
"  '  strands 

^  =  ^TtRN^trands  "  ^2^)^820)R  =  5152R  ft. 


873. 5R  KlVh 
Energ)''  Density  =  lb. 


watt -hr 
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Cost  =  .18$/ft.  X  5152R  ft.  =  $92 7R 

Cost  _  927R  $  _  $1062 

KlvTi  .8735RKlVh  KlVli 

Bearing  Power  Loss  (1-1/4"  shaft,  Fig.  9,  Ref.  4) 


E 


L 


.035 


watts 

lb-1000  rpm 


or 


watts  _  1  lb  _  ,  (-.j  v\fatts 

Ib-lOOOrpm  ^  .00773  KWli  KWh-lOOOrpm 


For  R  =  2 ' 


Store  1.747  Kl\fli  @  5480  rpm 

Weight  =  226  lb.;  10,304  ft  of  1/4"  rope 

Fl)n\,teel  Cost  =  $1855,  Bearing  Power  Loss  =  45  watts 


B.  KEVLAR^  hlie el  (1/4"  (|)Phillystran  Rope) 

For  8"  cross  section  (W  =  .021  Ib/ft) : 

A  =  ~  =  50.26  in^ 

^strand  ~ 

CJ  =  67501b/. 049  in^  =  137,755  psi 

\  =  (50 .26) ( . 8)  (density  factor) 

‘  strands  .049  ’  “  820  strands 

P  =  (820) (6750)  =  5.53  x  10^  lb. 


■R  = 


R’  I 


(2.654) (10^) 


=  (3.28)RKl\'h 


60 

27rR 


=  A  /  (2.768  X  10^)  (32.2) 
2rrR  V  (.021)  (820) 


21,735/R  rpm 
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W  =  2;rRN  ^  w  =  (2:7)  (820)  (.021  lb/£t)R  =  (108)R  lb. 
strands 

X  =  =  (2rT)(820)R  =  (5152)R  £t. 

3280R  (watt-h)  _ 

Energy  Density  =  ~(lb)~  "  lb 

Cost  =  X  5152R  £t.  =  3091R$ 


Cost  ^  30m  ^  0^1  3 

3.28R  KlVh 

Bearing  Power  Loss  (1-1/4”  shaft.  Fig.  9,  Ref.  4) 

£  =  Q3I^ - - 

\  lb- 1000  rpm 


r  _  n  ■’I-  watts  _1 _ 

\  -  U..^b  ^]3_x000rpm  ^  .03037  KlMi 


lb  ^  1  ^£7  watts 


KWh-lOOOrpm 


For  ^ 

Stores  6-1/2  Mi  @  10,862  rpm 

Weight  =  216  lb.;  10,304  ft  of  1/4”  rope 

Fl)nvheel  Cost  =  $6182,  Bearing  Power  Loss  =  82  watts 


C .  KElOJiR^  Wlieel  (1420  Den.  Kevlar-49)  F.S.  =  2 

3 

For  8"  cross  section  (density  (A)  is  .052  Ib/in  ): 


A  = 

P 


Tr8‘ 

~T 


Strand 


=  50.26  in^ 
=  68.81  lb 


68.81  lb. 

-4  .  2 
1.68  X  10  in 


409,600  psi 


N 


strands 


(50.26) (.9) 
(.000168) 


=  269,290  strands 


P  =  (269,290)  (68.81)  =  18.53  x  10  lb 


E  = 


'"-I 


N  = 

‘T  2^R 


(2.654) (10^) 
60 


=  (10.97)R  KlVh 


4)8 


w 


60 

2uR 


(18.53)(10  )(32.2) 
(2) (28.3) 


=  31s005/R  rpm 


=  (A) (A) (12) (.9)  =  28.29  Ib/ft 
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W  =  (x)CA)(2>rR)(12)(.9)  =  177. 35R  lb 


Encrov  Density  = 


10,970R  watt-h 
177, 55R  lb. 


,,  or  watt-h 

61.85  -yj— 


Cost  =  $8. 50/lb.  X  177. 55R  lb.  =  150 7R8 


Cost  _  150 7R 
KWIT  10.97R 


$157/KlVli 


Bearing  Power  Loss  (1-1/4"  shaft,  Fig.  9,  Ref.  4J 

^  _  orr  watts  ,,  1  =  |;7,1 

hi  Ib-lOOOrpm  .061  KlMi  '''  Klvh-lOOOrpm 

^  watts  „  1  lb  _  ^7,1  _ 

\  ~  Ib-lOOOipm  .061  KWh  KlMi-lOOOi’pm 


For  R  =  2' 

Stores  22  lOAli  @  15,500  I'pm 
Weight  355  lb. 

Cost  =  $3015  ('laterial  only) 
Bearing  Power  Loss  =  196  watts 
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APPENDIX  B  FLYlflEEL  EQUATIONS  OF  NOTION 


'Die  flwhecl's  di  fferent  i;il  ecjuations  of  motion  are  best  derived  by  way  of 
Layrangian  meduinics.  Ilie  general  form  of  Lagrange's  equation  of  motion  is 


_d  aL  .  aL  _ 

dt  ^  a^.  ’  a7[7  “  ^'i 


for  i  =  I ,  2 ,  ...  ,  n 


where  the  Lagrcrngictn,  L,  is  the  difference  between  the  kinetic  energy  of  the 
svstem,  T,  and  the  potential  energy  of  the  system,  V. 

L  =  T  -  V 

Ibe  is  the  generalized  nonconservative  force  acting  on  the  ith  coordinate,  ,  :md 
n  is  the  niunber  of  degrees  of  freedom. 

.■\n  exact  expression  for  the  system's  total  kinetic  energ\'  while  in  a  completely 
arbitran’  configuration  described  by  the  five  generalized  coordinates  sliovvn  in 
Figure  B1  is 


T  -  +  Llpu"  +  +  a“cos“S)  +  L(I^  +  sin'S 


C  A' 


+  LI^a)“  +  1^0001  sin  3  +  LI^((})“sin“(i)  +  cf)"')  +  LI ( (i)COSc|)  +  " 

where  Vp  is  the  translational  velocity  of  the  motor  mass  center,  G, . 

=  (h  +  OG^cosS)“a"  +  (OG^3)" 


,md  is  the  translational  velocity  of  the  flwheel  mass  center,  , 


2' 


Vp“  =  (h  +  OS  cos6)“a“  +  0S“3’“  +  fsin'sl"  +  £“0“  +  2£(h  +  OS  cos3)a 
(0cosesin((()-a)  +  ^sin  3  cos((l)-a))  +  2£  OS  3  sin  3  (4>sin0  sin((j)-a) 


-  0COS0COS  (l)-a) )  +  2£OSe6sin0cosB 


Tlrese  two  terms  in  T  express  the  overall  system  translational  kinetic  energy.  The 
remaining  terms  in  T  are  the  various  components  of  the  system  rotational  kinetic 
energy. 
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An  exact  expression  for  the  system's  total  potential  energy  is 
V  =  -Mggbcosa  -  (M^  +  M^^ghcosa  -  (M^OG^^  +  M2OG2)  geosgeosa 

2  2 

-  M2gfcos(j)sine  + 


The  various  length  parameters  appearing  in  these  expressions  are  defined  in  Fipre 
B1  except  for  b  which  is  the  distance  from  point  Q  to  the  center  of  mass  of  the 
horseshoe  shaped  gimbal  bracket.  The  various  properties  are  defined  as  follows: 


M-.  "  mass  of  the  motor,  including  any  shaft  and  any  motor  hub  to  which 
the  support  ropes  may  be  fastened 

-  mass  of  the  f  Indie  el 
NL  -  mass  of  gimbal  bracket 

D 

-  rotational  inertia  of  gimbal  bracket  about  point  Q 

I  “  rotational  inertia  of  the  motor's  non'rotating  parts,  mainly  the 
housing,  taken  about  the  motor's  main,  longitudinal  axis 

I  -  rotational  inertia  of  motor's  non-rotating  parts  about  a  transverse 
^  line  passing  through  the  motor  mass  center,  G2 

I,  -  rotational  inertia  of  motor's  rotating  parts,  mainly  the  armature, 

^  shaft  and  hub,  about  main  longitudinal  axis 

T  -  rotational  inertia  of  motor's  rotating  parts  about  a  transverse  line 
^  passing  through  G2 

I.,  -  rotational  inertia  of  flywheel  about  its  main  axis  of  rotation 

-  rotational  inertia  of  flindieel  about  a  transverse  line 
g  -  acceleration  due  to  gravity 


With  the  above  expression  for  the  Lagrangian  the  left-hand  side  of  Lagrange's 
equation  of  motion  can  now  be  obtained.  The  generalized  nonconservative  forces 
appearing  on  the  right-hand  side  are  obtained  from  an  expression  for  the  virtual  work 
done  bv  all  the  non conservative  forces  present  in  the  system  as  the  system  is  dis¬ 
placed^  through  a  virtual  displacement.  For  the  flywheel  system  the  non conservative 
forces  are  air  drag  acting  on  the  flywheel,  internal  hysteresis  in  the  support 
ropes  and  frictional  damping  on  the  two  gimbal  axes  of  rotation.  For  a  more  thoioug 
treatment  of  the  Lagrangian  method  of  analytical  mechanics  see  Reference  6. 
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A.  Rigid  Shaft  Model 


By  taking  the  derivatives  indicated  by  Lagrange's  equation,  along  with  the 
nonconservative  forces,  one  obtains  a  set  of  extremely  long  and  highly  non-linear 
differential  equations.  The  solution  to  these  equations,  if  it  could  be  obtained, 
would  be  an  exact  expression  for  every  possible  motion  the  system  can  have  under 
the  action  of  the  assumed  nonconservative  forces.  For  the  purpose  of  using  the 
flywheel  as  an  energy  storage  device  one  does  not  need  to  know  all  the  possible 
motions.  It  will  be  sufficient  to  know  whether  or  not  pure  spin  about  the  verti¬ 
cal  axis  is  a  dynamic  equilibrium  for  the  system,  and  if  so,  whether  or  not  this 
equilibrium  is  stable.  This  is  done  by  approximating  all  the  non-linear  terms  to 
be  linear  for  small  magnitudes  of  motion  about  the  equilibrium  configuration  in 
question.  This  will  reduce  the  set  of  non-linear  differential  equations  to  a  set 
of  linear  differential  equations  whidi  can  be  solved  by  applying  linear  oscillation 
theory . 

For  the  rigid  shaft  model  this  set  of  linearized  differential  equations  for 
the  case  of  constant  motor  speed,  o,  is 

I^a  +  (h  +  0S)$  +  I.og  +  K  a  =  (h  +  OS) (F^  +  F^)  -  C  a 

z  A  a  a  a  a 

IqB  +  ose  -  +  Kg3  =  -OS  (F^  +  F^)  -  C 

+  ILl  (h  +  0S)d  +  I„w0  +  K  (j)  =  £F^ 

1  Z  Z  (p  (f) 

i5e  +  M.,  ose  -  1 70)0  +  K„0  =  -tf!’ 

1  Z  Z  0  0 


Each  term  in  these  equations  represents  a  torque  that  is  present  in  the 
system  such  as  inertia  torques,  gyroscopic  torques,  gravity  torques  and  non¬ 
conservative  torques  due  to  friction.  The  constants  appearing  are  length  para¬ 
meters  shown  in  Figure  B1  and  system  inertias  and  springs  defined  as  follows: 

Iq  =  +  l£  ,  effective  system  rotational  inertia  about  upper  gimbal 

axis  through  point  Q 

0  —  2  —  — 

m  ~  ^  ^B’  rotational  inertia  of  motor  and  support 

system  about  upper  gimbal  axis  through  point  Q 

Q  _  2 

r  ~  ^  ’  effective  rotational  inertia  of  flywheel  about  upper 

gimbal  axis 

Iq  =  +  l£  ,  effective  system  rotational  inertia  about  lower  gimbal 

axis  through  point  0. 

0  _  2  —  — 

,  rotational  inertia  of  motor  and  support  system 
about  lower  gimbal  axis. 
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=  M^OS^,  effective  rotational  inertia  of  flywheel  about  loiver  gimbal 
axis 
S  2 

^f  ~  ^  inertia  of  flywheel  about  a  horizontal 

line  through  point  S 

=  (Mgb  +  M^h  +  +  M^OG^  +  M20S)g  +  ,  total  spring  rate  asso¬ 

ciated  with  the  a  coordinate 

k  ,  spring  rate  of  any  spring  that  is  installed  in  addition  to  gravity 
“  on  the  a  coordinate 

Kg  =  (M^OG^  +  M20S)g  +  kg  ,  total  spring  rate  associated  with  the  g 
coordinate 

k  ,  spring  rate  of  any  spring  installed  in  addition  to  gravity  on  the 
^  g  coordinate 

gra\/ity  spring  rate  on  tj)  coordinate 
Kg  =  M-^fg,  gravity  spring  rate  on  9  coordinate 

Tlie  components  of  tlie  nonconservative  forces  are  defined  as  follows: 

C  ,  viscous  damping  constant  for  motion  in  the  a  coordinate  (torque 
™  per  angular  velocity) 

C  ,  viscous  damping  constant  for  motion  in  the  g  coordinate  (torque 
^  per  angular  velocity) 

=  -C^[co(-£e  +  £g)  -  {t\  -  fa)],  internal  friction  force  acting 
at  the  motor  hub,  point  S,  in  the  a  coordinate 

C.  ,  internal  friction  constant  of  proportionality  (force  per  linear 
velocity) 

0  =  e  -  I 

Fn  =  =  -C.[-w(f(|)  -  fa)  -  (-f0  +  fg)],  internal  friction  force  acting 

g  g  1 

at  the  motor  hub ,  point  S ,  in  the  g  coordinate 

pG  _  ipG  ^  apG  -total  nonconservative  force  acting  at  fl wheel  center, 

(f)  (})  (}) 

G^^  in  (f)  coordinate 
i  G  .  . 

F  =  -Gitif.  (4)  -  a)  +  C.f(6  -  g) ,  internal  friction  force  at  flywheel 
0  1 

center  in  coordinate 

Q  ,  _ , 

F  =  C  (fe  +  OSg) ,  air  drag  force  at  fl>^vheel  center  in  6  coordinate 

0  3. 
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The  difference  between  the  bowed-out  rope  model  and  the  rigid-shaft 
model  is  the  configuration  of  the  support  ropes.  This  difference  causes  two 
basic  changes  in  the  differential  equations  listed  above.  First  the  ilywheel 
will  now  behave  as  if  it  were  a  point  mass.  Second,  the  different  geomet^ 
of  the  support  ropes  changes  the  action  of  the  internal  friction  force.  These 
changes  are  summarized  as  follows: 


=  -C.aj£(-e  -  fg)  +  C.lC^  -  fa) 
a  1  1 

^Fg  =  C^a)f,(<|)  +  fa)  -  C^O)  (0  +  fg) 

^F^  =  C.(jil<i-e  -  fa)  -  C.  (4)  -  fa) 
a  1  1 

^Fg  =  -C^caf((t)  +  fa)  +  (6  +  fg) 

where  the  constant  f  is  defined  as 
f  =  a/(D  -  a) 

with 

a  =  diameter  of  the  motor  hub  where  the  support  ropes  are  attached 

p  =  diameter  of  the  flywheel  where  the  support  ropes  are  attached 

With  these  changes  noted  the  four  differential  equations  of  motion  for  the 
bowed-out  model  are  identical  to  those  derived  for  the  rigid-shaft  model. 
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APPENDIX  C 


Computer  Program  for  Stability  Analysis 


The  eigenvalues  and  eigenvectors  for  the  matrix  A  of  equation  (1)  are 
found  by  using  an  International  Mathematical  and  Statistical  Library  (BfiL) 
program  by  the  name  EIGRF,  The  program  is  designed  to  find  the  complete 
eigenvalues  and  (optionally)  eigenvectors  of  a  real  general  matrix  in  full 
storage  mode.  The  algorithm  used  by  this  program,  see  References  (1,2) 
makes  use  of  a  number  of  other  IMSL  subroutines.  The  matrix  is  first  bal¬ 
anced  to  reduce  numerical  error  and  then  a  Hessenberg  matrix  is  found  for 
this  balanced  matrix.  The  eigenvalues  and  eigenvectors  are  found  for  the 
Hessenberg  matrix  and  are  then  back -transformed  into  the  eigenvalues  and 
eigenvectors  for  the  original  matrix  A. 

Once  the  eigenvalues,  X,  and  the  eigenvectors,  u,  have  been  found,  the 
stability  of  the  system  can  be  deduced  directly  from  the  eigenvalues. 
However,  to  find  the  mode  shapes  from  the  conplex  eigenvectors  requires 
some  manipulation.  First,  it  should  be  noted  that  since  matrix  A  contains 
only  real  numbers  it  can  be  shown  that  the  eigenvalues  will  always  appear  in 
complex  conjugate  pairs,  and  from  this  it  can  be  shown  that  the  correspond¬ 
ing  eigenvectors  will  also  be  complex  conjugates.  Knowing  this,  the  general 
solution  for  the  state  vector,  X,  can  be  written  as 

X  =  B,u, e^^  +  B^u,e^l^  +  . . .  +  B^  ^u  e^n^  +  B^  u  e^n^ 

11  21  2n-l  n  2n  n 

where  the  B's  are  complex  constants  and  the  bar  denotes  complex  conjugate. 

One  complex  conjugate  pair  of  eigenvalues  contributes  two  independent 
parts  to  the  general  solution,  which  is  also  complex.  But  in  the  case  of 
the  flywheel,  which  is  a  real  system,  the  solution  needs  to  be  real  so  that 
we  don't  want  the  entire  general  solution.  For  the  solution  to  be  real, 
the  constants  associated  with  a  conjugate  pair  of  eigenvalues,  say  B,  and 
B^,  must  then  also  be  complex  conjugates.  This  will  make  these  two  terms 
of  the  general  solution  complex  conjugates  of  each  other  and  therefore  yield 
a  real  expression  upon  addition. 

There  are  still  two  arbitrary  constants  associated  with  the  pair  of 
eigenvalues,  the  real  part  and  the  imaginary  part.  To  find  the  mode  shape 
associated  with  a  pair  of  conjugate  eigenvalues ,  set  the  real  part  of  the 
constant  equal  to  one  and  the  imaginary  part  equal  to  zero.  This  will  give 
real  numbers  for  the  coordinates  of  the  system,  and  if  desired  can  be 
plotted  to  show  the  relative  whirl  orbits  of  the  flywheel  and  the  motor 
which  in  general  will  be  elliptical. 

A  listing  of  the  calling  program  follows: 
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(JUWf 
T  T  C  ^  ; 

FT  CT  : 

^  '  0 

34 

3  0 


7P)  OS/SfcO  FORTRflKHrTFKDED  rAfE79o 

^  no  F  C  F  9  N  (i  L  1ST  9  C  P  T  (  0  J  9  K  0  D  U  F  P 

NftMr(MAlN,)  KiOOP  T  1  LI  ME  C'^U  ^T^  6  0  )  SIZEd'^AXli  f- U  t  f  D  P  L  «  M  0  N  E  ) 

SOURCE  rOCDIC  UCLIST  KODECk"  OE.JECT  NOF'AP  AfFCRRAT  GCSTr^T  \O^REF  AlC  NOAE'S 
(;  o  M  D  I  r  i  5  ^  L  A  W  P  n  A  «  P  )  9  2  (  <5  H  )  9  Z  N 

E  0  U  'U.  r.  PR  E  C  I  S  I  C  El  E  <  P  9  8  1  9M  1  •  R?  9  0  9  1  P  9  I  A  9  I  C  A  9  I  ■  9  I  2  9  L  9  H  <3  I  HG  9  1  F  E  9  I  S  9  I 

l!O9lf'':9l0  9l'*S9KA9KB9t<'PHloKTEE9HCG9H('S9BR»0G90S9O.  RAVQCA9CB9CI9CO9CID 
‘t  9  w  9  A  9  A  r  9  A  P  P  ,0-9  P  P  9  B  F  P  9  C  9  CP  9  C  PP  9  r  9  DP  9  DP  P  9  E  K  C  1  P  0  )  9  r 
data  e  /  (,  4  *  c  o  r  0  / 

P  r  A  P  (  S  9  3  0  0  )  1^  0  1  9  ?  9  1  R  9  I  A  9  I C  A  9  r  9  I  2  9  H  9  P  P  9  0  F  9  0  S  9  L  9  G  R  A  V  9  C  A  9  C  E'  9  c  I  9  C  C  9 

J?  'a' 

Frii5‘^AT{3FlGoE9/9  5Flt'o59/9feFlOo'=''9/94F10oG9/9riEo'^5 
DO  '3  I  I  =  1  9  £ 

ri  =  E.PP4D090<,0  H'lO^JFLOATJl  1)-1,0) 

WRITr«t9’0D)MpoMj,wp,lr.  ,lA,TCA9lA'9l29H9-R,CG9CS9L9GRAV9CA9Cri9Cl9CD 

q 

W  J  T  r-  (  e  9  3  6  0  ) 

r  ^  f  1  9  /  /  /  /  ) 

H  0  G  r  H  r  G 
l-J  n  G  "  G 

(HOG^HOG^^-^ICA^IF 

1  r  -  M  2  ^  f  H  0  0  > 

2  Q~  l^Q  Q 
IPG-F^l  *COG*CG!>  <^ICA 

1  r  0  -  2  *  0  s  0  s ) 

TO-'l'^o^iro 

ypCr»^2trL^L’^T 

KA:r(M.  *OG<^i^G«OG)'^iGPAV 

KP-  ?  ^og^i^g^ogmgpav 
KP'HJ  tM  2  A  gravel 
K  THG  rhi 

jr  j  i  \  HOS^  tCS  y 

AP-:rS/fM2*L^HGS5 

BP-irs/<M2‘'L*o:^) 

C-- A 

CP  -  T  G  /  (  1^2  *  L  *  HO  *': ) 

D  -  R 

DP-If^/^H2*L*CS> 

r=io87cnc/«poOrc-ioe.75r:G5 

CCl<5l)~<AP^CD^^hOS^HGG~(:D*l  ^HOS“CI^L’^^L‘«=^r^AP»aCA  )/A 
rOgD-AP^rjSt^W^A 

E  (  1  9  3  >  ^  ?  AP  ^HOS  *  C  D  ^  C  I  <>CD  )  *L  L  )  /  A 
r  (  1  9  4  )  -  C  o  C  c: 

ECloG)-AP*PA/A 
':(ioE)-“«ci^w^L^L*r/A 
rCl975~-KPHT/A 
r  i  ioR>=-rT*W^L^L/A 
E  ?  r  0  1  )  =-RP  *  T  A  *  G/R 

E<2QAG>  =  CBP’^Cr^rG<^0E-CD*Lt^C£”>CI^L*L^r«pr-^CB>/B 
E  (  2  9  3  >  -  0  o  0  D  0 

r(2^4)  r(-^CT^Cr>^^L’!^L^8c^*0S  H.  C  C  >  /B 
E  ?  P  9  5  )  =  C  I  W  ^  L  L  r  /  B 
E(P<5G>-BP*^G/B 
E  C  P  0  T  }  -  C  I  ^  W  L,  ^  L  /  B 
r  ^  2  qR  ) -«K  THE /R 

r  (?Q  1)  ~(-rofrMOE*Hos^cp*ci  L’^l^F'S^cp^cd^l^^hcg-ca  j/c 
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(JUNt  7(") 


INI 


OS/ 36  0  FORTRAN  H  F 'TENDED 


DATE  79 


E  <  ^  ,  2 ) =-I A  *  W/C 

r<3*'5)  =  (-  CD»l*HOS-»CP<*L*L*(CI*CD))/C 
E  (  3  0  A  )  =  0  .  C  D  0 
E ( 3 ,E ) --« A /C 
E  (3»6)  =  CP*ri*W*L*L»r/C 
»  E(3*7)=CP*<PHI/C 

E (3«P>-EP*CI*W  *L*L/F 
F  (  A  ^  j ) -lA  *  w/D 

=  CD*0S»0S+DP*CI*L*L*F*nP*CD«L*0‘^ 
’  C(A,2>rO.OCn  u. 

E (A  »4  )  =  (DP* (Cl  +CD) *L*L-CD*  L*OS  I/D 
E  (A9E)--DP»CT*  i«*L*L*F/C 
E  (  A  ,  <•  )  --KP  /D 
E (A  V  7)  =-DP  *CI*W*L*L/D 
E(4,“)=DP*KTHE/D 

E  (  6  .  r? )  - 1 . 0  D  0 
E(7,3>=1.0DC 
r(P*A):r3,0D'' 

Ki=U 

IJ  0  9  =  2 

E-ALL  EIGRF  (E»N?N*lJOB*LAMBEAisiZ  ?N‘»W(<«irR) 
DO  E  J  =  1,N 
3  N>  -  2  (  7  <  O  ) 

DO  9  I  :r  1  ,  N 
2 {  I  ,  J) -2(  I  ,J  )/2N 
E  CONTirjLE 

DO  6  I =  3  «  M 

WP  I  TE  {  t  ,1  OC  )  LA  NBDA  (  I)  ,  (C(  I  «vJj  ,  Jrrl  ,  N  ) 

I'O  FORVAT(1a»fio»?:»3v,di?.i2,»  a9,FDi3.E) 

6  0  0  N!  T  I  I.  £ 

WPlTr(6,12E) 

125  FORPAT  <1  :■  ,////  ) 

DO  7  I =  1  ,  N 

UiPTTt  (6,150)  (Z  (I,  J)  ,J=1,4  ) 
lEO  FOR y A T ( 1 V , poi 3 *5 > 

C  0  V  T  I  M  I,  F 
WRITE  (6,175) 

175  FORf'ATflv,//) 

DO  R  I  =  3 ,  P 

WRITr-(fc,i50)  (Z  (I,  J)  ,J='-',fi) 

6  CON'TIN'Ur 

WRlTr{6,125) 

URITE(6,200)  WK  (1),  lER 

2(  0  FOPNAT  (1 V,////,D2  3.16,5  ,  I  E,/, IHl ) 

13  C  C  '\l  T  I  NJ  L  r 

STEP 
F  N  D 


-CP)/D 
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